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Abstract 


Objectives 

An  essential  eomponent  in  the  deteetion  and  charaeterization  of  underwater  munitions  is 
knowledge  of  the  acoustic  response  of  the  environment  as  well  as  the  environment's 
effect  on  the  acoustic  response  of  munitions.  Simulation  tools  and  technologies  have  been 
developed  under  SERDP-funded  research  initiatives,  such  as  the  Personal  Computer 
Shallow  Water  Acoustic  Toolset  (PC  SWAT),  to  model  the  acoustics  of  both  the 
environment  and  the  munitions.  The  evaluation  of  these  technologies,  as  well  as  their 
future  use  in  unexploded  ordnance  (UXO)  remediation,  rely  on  knowledge  of  the 
underwater  environment,  particularly  the  properties  of  the  seafloor.  Conventional 
methods  for  determining  relevant  seabed  properties  employ  time-consuming  point 
measurements.  Commercially  available  high-frequency  multi-beam  echo  sounders 
(MBES)  offer  a  solution  to  this  problem  by  making  faster  measurements  over  larger 
areas.  While  these  systems  are  primarily  intended  for  high-resolution  bathymetry, 
acoustic  inversion  techniques  can  be  used  to  estimate  seafloor  parameters  relevant  to  the 
UXO  problem.  The  objective  of  this  research  is  to  develop  and  rigorously  test  a  physics- 
based  algorithm  that  can  invert  high-frequency  acoustic  data  for  sediment  parameters. 
The  result  is  a  collection  of  new  algorithms  for  high  frequency  acoustic  data  inversions 
and  system-independent  environmental  assessment  in  terms  of  measurable  seabed 
parameters,  which  can  then  be  used  as  inputs  to  acoustic  and  electromagnetic  systems. 

Technical  Approach 

The  inversion  algorithms  developed  here  incorporate  the  results  of  recent  research  in 
high-frequency  sediment  acoustics,  utilizing  both  models  and  technologies  that  have  been 
developed  through  ONR-funded  experiments  such  as  the  Sediment  Acoustics 
Experiments  1999  (SAX99)  and  2004  (SAX04).  To  develop  and  test  the  algorithms, 
acoustic  data  was  collected  during  three  field  experiments: 

1 .  The  Gulf  of  Mexico  Experiment  in  201 1  (GulfExl  1), 

2.  The  Target  and  Reverberation  Experiment  in  2013  (TREX13), 

3.  The  St.  Andrew’s  Bay  Experiment  in  2014  (BayExl4). 

These  field  efforts  were  sponsored  by  both  SERDP  and  ONR  under  several  different 
programs  and  this  project  utilized  the  infrastructure  and  logistics  of  those  efforts  to 
collect  data.  In  each  of  these  efforts,  the  acoustic  data  was  collected  using  Seabat  7125 
multibeam  echo  sounders  that  were  provided  along  with  technical  support  by  Teledyne- 
RESON.  The  acoustic  data  collection  in  each  of  these  experiments  was  accompanied  by 
extensive  environmental  characterization  in  order  to  evaluate  the  accuracy  of  the 
inversion  algorithms. 
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Results 

While  the  data  collected  during  GulfExl  1  provided  the  basis  for  developing  the  inversion 
algorithm  and  to  learn  what  was  needed  to  collect  a  quality  dataset,  equipment  problems, 
the  abundance  of  fish,  and  saturation  of  the  MBES  data  prevented  the  use  of  that  dataset 
for  inversion  testing.  The  lessons  learned  during  this  experiment  were  applied  to 
TREX13  and  BayExl4  and  helped  to  insure  that  the  data  was  of  sufficient  quality  to 
evaluate  the  inversion.  The  environments  in  both  of  these  experiments  were  complicated, 
but  the  inversion  performed  well  producing  stable  results  that  could  be  compared  to 
ground  truth  measurements  at  the  sites.  These  ground  truth  measurements  were  either 
compared  directly  to  the  inversion  outputs  or  used  in  scattering  models  to  compare  to  the 
inversion. 

For  TREX13,  the  inversion  compared  favorably  to  the  ground  truth  measurements  for  a 
majority  of  the  data  products.  Discrepancies  between  the  scattering  models  and  the 
inversion  are  most  likely  due  to  the  inability  of  the  scattering  models  to  account  for  the 
multiple  scattering  from  shell  pieces  at  the  site.  The  roughness  spectrum  obtained  from 
the  inversion  compared  well  with  the  measurements  across  the  entire  test  area  for  high 
wavenumbers  indicating  that  the  results  would  be  applicable  to  modeling  shallow  grazing 
angle  scattering  for  frequencies  at  least  as  low  as  50  kHz. 

The  environment  at  the  BayExl4  site  was  significantly  more  complicated  than  the 
TREX13  site.  The  presence  of  the  mud  layer  and  the  transition  layer  between  the  mud 
and  sand  sediments  was  not  accounted  for  in  the  inversion,  which  treated  the  seafloor  as  a 
half-space.  Despite  this  approximation,  the  inversion  results  are  stable  and  give  effective 
parameters  for  the  sediment  that  differ  from  the  ground  truth  measurements.  The 
scattering  models,  which  incorporate  the  ground  truth  measurements,  also  do  not  account 
for  the  complexity  of  the  environment  and  could  potentially  be  modified  to  more 
accurately  represent  the  scattering  physics. 


Benefits 

Over  the  course  of  this  project,  the  MBES  inversion  has  reached  a  level  of  maturity  that 
puts  it  at  the  cutting  edge  of  remote  seabed  characterization  particularly  for  the  EIXO 
detection  problem.  Taking  a  physics-based  approach  insures  that  the  output  of  the 
inversion  addresses  the  needs  of  modeling  codes  and  simulations  that  incorporate  similar 
physical  models.  This  inversion  algorithm  and  the  data  collection  approach  needed  to 
apply  it  have  the  potential  of  being  implemented  into  sonars  without  modifying  existing 
hardware  or  software. 
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Objective 

The  objective  of  this  research  is  to  develop  a  technique  that  can  be  used  with  multibeam 
sonars  to  provide  rapid,  wide-area  assessment  of  seafloor  properties  that  can  aid  in  the 
detection  and  classification  of  UXOs.  While  electromagnetic  techniques  have  shown 
strong  promise  as  a  means  of  detecting  and  discriminating  UXO  from  rocks  and  other 
debris,  they  must  be  close  to  the  target  in  order  to  work  well.  This  limits  their  capability 
in  wide-area  surveys  of  remediation  sites  and  will  need  to  be  augmented  by  other 
technologies.  Acoustic  systems  such  as  synthetic  aperture  sonar  (SAS)  may  provide  a 
complementary  means  of  detecting  UXO  over  a  wide  area,  but  the  ability  of  this 
technique  to  discriminate  and  classify  UXO  is  heavily  influenced  by  the  underwater 
environment  in  which  they  are  operating.  Sediment  roughness  and  volume  heterogeneity 
both  contribute  significantly  to  bottom  reverberation,  the  dominant  limiting  factor  for 
UXO  detection  and  classification.  Reverberation  appears  as  a  background  in  the  total 
acoustic  response,  from  which  the  target  contribution  is  to  be  detected  and  classified.  The 
ability  of  the  sonar  propagation  and  scattering  models  such  as  the  PC  SWAT  to 
accurately  model  the  reverberation  depends  on  knowledge  of  the  sediment  properties. 

Conventional  techniques  to  measure  sediment  properties,  such  as  stereo  photography  and 
sediment  sampling,  are  time-consuming  and  labor-intensive.  They  also  provide  only 
point  measurements,  and  multiple  measurements  are  required  to  provide  sufficient 
information  about  a  complex  environment.  Remote  sensing  techniques  using  sonar  have 
long  been  seen  as  a  means  of  providing  information  about  the  seabed  over  large  areas  at 
potentially  high  resolution.  The  drawback  to  most  implementations  of  sonar  for  seabed 
characterization  is  that  they  still  rely  on  ground  truth  measurements  to  “calibrate”  the 
algorithms  used  to  identify  seafloor  types.  In  the  end,  these  techniques  don’t  replace  the 
point  measurements,  but  instead  supplement  them  by  taking  the  information  gleaned  from 
the  point  measurements  and  mapping  it  over  a  wider  area. 

To  truly  be  time  and  cost  effective,  the  acoustic  seabed  characterization  should  work 
independently  from  any  ground  truth  measurements.  Evaluation  of  the  sediment 
properties  would  then  involve  a  single  acoustic  survey.  For  this  approach  to  be  viable, 
the  acoustic  inversion  would  need  to  be  system  independent  and  physics-based.  Research 
into  high-frequency  acoustic  interactions  with  the  seafloor  have  progressed  significantly 
over  the  last  30  years  and  has  reached  a  level  of  maturity  where  this  type  of  physics- 
based  approach  may  be  feasible.  The  objective  of  the  research  reported  here  is  to  develop 
such  a  physics-based  algorithm  and  evaluate  its  performance  using  MBES  data  collected 
at  three  experiment  sites.  The  ultimate  goal  is  to  develop  the  algorithm  to  a  point  where 
it  can  be  integrated  into  the  sonar  processing  suites,  extending  the  capability  of  existing 
bathymetric  sonars  to  include  seafloor  characterization.  This  type  of  system  could 
address  not  only  SERDP’s  needs  but  also  provide  a  valuable  tool  for  marine  habitat 
mapping  and  hydrographic  surveys. 


3 


Distribution  Statement  A:  Approved  for  public  release;  distribution  unlimited. 


Background 

Traditionally,  the  nature  of  the  seafloor  has  been  determined  by  means  of  direet 
sampling,  photography,  video,  and  various  in-situ  probing  methods.  Sueh  measurements 
are  time  eonsuming  and  expensive  when  a  large  area  is  to  be  eharaeterized.  Aeoustie 
methods  have  been  developed  with  the  aim  of  providing  rapid  eoverage  of  large  areas. 
Two  eommereially  available  methods  have  been  widely  used  for  over  a  decade.  Though 
variants  are  available,  these  approaches  are  exemplified  by  QTCViewTM™  [l]-[3]  and 
Rox  AnnTM™  [4]-[6].  These  methods  extract  “features”  from  the  sonar  echo  envelope 
and  subject  these  to  clustering  analysis  to  classify  portions  of  the  survey  area  into  several 
types.  At  this  stage,  the  true  nature  of  these  types  is  unknown,  as  sonar  echoes  depend 
not  only  on  the  properties  of  the  seafloor,  but  also  on  water  depth  and  sonar  properties. 
In  order  to  identify  clusters  with  specific  seafloor  types,  acoustic  surveys  must  be 
accompanied  by  ground  truthing  using  methods  such  as  grab  sampling  and  video 
observation.  If  sufficient  ground  truth  is  available  to  include  all  seafloor  types  in  the 
survey  area,  the  clusters  obtained  from  acoustic  data  can  be  identified  with  seafloor  types. 
Where  time  and  resources  have  permitted  sufficient  ground  truthing,  these  methods  have 
proven  to  be  robust  and  useful. 

As  ground  truthing  requires  extra  resources  and  can  be  expensive  and  time  consuming, 
there  have  been  efforts  to  avoid  this  step  by  extracting  echo  features  that  are  intrinsic  to 
the  seafloor  and  independent  of  such  factors  as  water  depth  and  sonar  parameters.  These 
methods  employ  physical  models  for  acoustic  scattering  by  the  seafloor  [7]-[10].  The 
power  of  this  approach  is  increased  when  used  with  multibeam  sonars,  which  typically 
provide  more  than  100  time  series,  one  for  each  beam,  greatly  reducing  inversion 
ambiguity. 

The  primary  challenge  in  the  approach  presented  here  is  that  the  models  used  must  be 
sufficiently  realistic  and  general  to  emulate  a  wide  range  of  seafloor  types.  Our  use  of  this 
approach  in  the  SERDP  application  is  supported  by  the  successful  development  of  a  wide 
range  of  seafloor  acoustic  models  by  personnel  of  APL-UW  [11]. 
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Materials  and  Methods 

In  order  to  develop,  test,  and  evaluate  the  sonar  inversion  algorithm,  we  collected  both 
acoustic  data  for  inversion  and  ground  truth  data  to  test  the  performance  of  the  inversion. 
The  systems  used  to  make  these  measurements  as  well  as  the  inversion  algorithm  itself 
are  described  below. 

Teledyne-RESON  Seabat  7125 

The  physics-based  inversion 
algorithm  developed  during  this 
project  was  not  designed  for  use 
with  a  specific  sonar  system.  It 
instead  requires  that  the  sonar  be 
well  characterized  (calibrated 
source  level,  beam  patterns,  and 
well-understood  pulse 

characteristics)  and,  for  accurate 
results,  have  narrow  beams  and 
short  pulse  duration  such  that  the 
interface  and  volume  response  of 
the  sediment  can  be  well  separated. 

In  order  to  collect  data  to  test  the 
inversion  algorithm,  for  this  project 
we  chose  to  focus  on  the  Teledyne- 
RESON  Seabat  7125.  This  was  a  sonar  system  that  has  been  used  in  other  projects  at 
APL-UW  and  by  working  with  our  colleague  Dr.  Gorm  Wendelboe  at  Teledyne-RESON, 
we  were  able  to  use  R&D  versions  of  the  7125  for  each  of  the  field  tests  at  no  cost  to  the 
project.  These  systems  had  modifications  and  improvements  that  made  them  ideally 
suited  for  the  data  needed  for  this  project  and  these  improvements  to  the  Seabat  7125 
have  since  been  incorporated  into  production  models  of  the  sonar. 

The  Seabat  7125  consists  of  a  source  and  receive  array  as  shown  in  Figure  1.  The  source 
projects  a  beam  pattern  on  the  seafloor  that  is  very  narrow  along  track  and  very  broad 
across  track.  The  receive  array  is  oriented  perpendicular  to  the  source  and  the  receive 
signal  can  be  beam-formed  to  produce  256  beams  each  of  which  is  broad  along  track  and 
narrow  across  track.  As  a  result,  the  signal  measured  for  a  given  beam  corresponds  to  the 
acoustic  response  of  a  small  portion  of  the  seafloor.  Typically  the  time-of-arrival  for  that 
response  is  used  to  compute  the  distance  to  that  patch  of  the  seafloor  and  then  used  to 
produce  high-resolution  bathymetric  maps.  The  full  time  series  response  of  each  ping 
can  also  be  recorded  as  was  done  in  each  of  the  experiments  discussed  here  and  used  in 
the  inversion  algorithms. 

The  Seabat  7125  is  usually  operated  at  either  200  kHz  or  400  kHz  depending  on  the  depth 
of  the  seafloor.  The  beam  width  at  200  kHz  is  2  degrees  along  track  and  1  degree  across 
track.  The  beam  width  at  400  kHz  is  1  degree  along  track  and  0.5  degrees  across  track. 
The  beams  extend  across  a  full  swath-width  of  140  degrees.  For  this  project,  the  software 
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Figure  2:  The  IMP2  (left  image)  and  the  SLS  (right  image), 

used  to  run  the  7125s  was  modified  to  operate  at  20-25  kHz  inerements  from  180  kHz  to 
420  kHz.  The  beam  widths  at  each  frequency  varied  between  the  two  extremes  at  200 
and  400  kHz. 

Environmental  Characterization 

While  the  inversion  algorithm  does  not  rely  on  ground  truth  data  to  identify  the  seafloor 
types,  evaluation  of  the  algorithm  performance  does  require  a  separate  measurement  of 
the  seafloor  properties.  This  seafloor  characterization  was  conducted  during  each  field 
test  to  try  using  systems  and  techniques  developed  over  30  years  of  high-frequency 
acoustics  research  by  APL-UW  and  our  collaborators. 

In-situ  Measurement  of  Porosity  2  (1MP2) 

The  IMP2  is  a  sediment  conductivity  measurement  system  developed  to  support  SAX99 
and  SAX04  [12],  [13]  (Figure  2).  It  uses  a  single  constant-current  electrode  that  is 
scanned  through  the  sediment  in  a  2D  plane  using  a  horizontal  and  vertical  positioning 
system  mounted  on  a  frame  that  rests  on  the  seafloor.  The  probe  can  be  inserted  up  to  25 
cm  into  the  sediment  in  1  mm  increments  and  at  1  cm  horizontal  steps  up  to  4  m.  At  each 
position  the  output  voltage  of  the  probe  is  recorded  which  is  inversely  proportional  to  the 
mean  conductivity  over  the  resolution  cell  of  the  electrode.  The  “formation  factor,”  F,  is 
found  by  taking  the  ratio  of  the  voltage  measured  in  the  sediment  to  the  voltage  measured 
in  the  overlying  water.  From  the  formation  factor,  the  sediment  porosity  is  calculated 
according  to  an  empirical  relation  known  as  Archie’s  Law, 

.  ^  J_  (1) 

^  pl/m’ 

where  m  is  a  constant  for  a  given  sediment.  This  constant  is  determined  by  calibrating 
the  probe  using  a  sediment  sample  that  is  later  dried  and  weighed  to  determine  the  bulk 
porosity.  The  spatial  distribution  of  the  porosity  fluctuations  can  then  be  calculated  and, 
if  the  sediment  material  properties  are  known,  the  density  fluctuations  can  also  be 
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calculated.  These  are  used  to  determine  the  fluctuation  power  speetrum,  whieh  is  used  as 
an  input  to  perturbation  theory  for  volume  seattering  from  the  sediment. 

Seafloor  Laser  Scanner  (SLS) 

The  SLS  is  a  struetural-light-based  system.  It  consists  of  a  waterproof  laser  line  projeetor, 
charge-eoupled  device  (CCD)  camera,  and  PC  eontrol  unit  [14]  (Figure  2).  The  scanning 
head  assembly  maintains  the  relative  position  and  orientation  between  the  eamera  and  the 
laser  line  souree.  The  eamera  is  a  Easier  A102f  CCD  with  a  resolution  1388  x  1038 
pixels.  The  laser  wavelength  is  635  nm  and  a  cylindrieal  lens  is  plaeed  at  its  tip  to 
generate  a  60-degree  fan  angle. 

The  eamera,  projeetor,  and  housing  are  mounted  on  the  IMP2.  The  SLS  scanning  head 
assembly  is  mounted  on  the  horizontal  linear  stage,  roughly  75  em  above  the  seafloor,  to 
achieve  an  effeetive  seanning  swath  width  of  30  em.  The  optical  resolution  on  the 
scanning  plane  (across  track)  is  0.3  mm.  In  each  image  frame,  the  high-contrast  pixels 
eorrespond  to  the  reflection  of  the  laser  line  from  the  seafloor.  With  a  calibrated  camera, 
the  pixels  on  the  laser  line  can  be  converted  into  the  aetual  dimension  of  the  bottom 
profile  at  the  scan  location  [15].  As  the  IMP2  linear  stage  moves,  the  laser  is  scanned 
along  the  traek.  The  resolution  along  the  traek  is  approximately  0.5  mm.  This  results  in  a 
digital  elevation  map  (DEM)  with  an  area  of  about  30  em  x  360  em  with  0.3  and  0.5  mm 
resolutions  across  and  along  the  traek,  respeetively,  and  0.1  mm  in  the  vertieal  direetion. 

The  DEM  are  used  to  ealeulate  the  ID  roughness  power  speetra  from  which  the  3D 
power  speetra  ean  be  determined  by  assuming  that  the  roughness  is  isotropie.  These  3D 
power  speetra  are  then  used  as  inputs  for  the  small-slope  interfaee  scattering  model  for 
comparisons  with  the  inverted  interface  scattering  strengths. 

Sediment  Sampling 

During  each  experiment,  both  diver  eores  and  diver  eofferdam  samples  were  colleeted 
and  analyzed.  The  diver  cores  are  5. 9-em- inner-diameter,  polycarbonate  plastic, 
cylindrical  cores  that  are  manually  pushed  into  the  sediment  by  a  scuba  diver.  The  top 
and  bottom  of  the  cores  are  eapped  and  brought  to  the  surface  for  later  analysis.  These 
analyses  inelude  measurement  of  the  porosity,  sediment  density,  and  shell  size 
distribution.  Porosity  was  found  by  measuring  the  water-weight  loss  of  the  sediment  by 
either  drying  the  whole  eore  or  a  seetion  of  the  eore  at  105  degrees  C  for  24  hours.  The 
sediment  density  was  similarly  analyzed  by  extraeting  the  excess  water  from  eore, 
weighing,  and  accounting  for  the  weight  of  the  eore  and  caps.  The  shell  size  distribution 
was  assessed  as  a  function  of  depth  by  sectioning  cores  and  dry  sieving  the  material  in 
eaeh  seetion.  The  sieve  sizes  covered  the  range  1  mm  to  9.5  mm.  The  material  at  a  eaeh 
sieve  size  was  then  weighed  and  eounted  to  determine  both  the  shell  weight  distribution 
and  the  number  density. 

Sinee  the  eores  eould  only  capture  shell  sizes  less  than  5.9  em  in  size,  larger  cofferdams 
were  also  used.  The  cofferdams  eonsisted  of  a  plastie  eylindrical  pipe  with  a  21-cm- 
inner-diameter  that  eould  be  inserted  10  em  into  the  sediment.  The  material  within  the 
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Figure  3:  The  "attenuation  array"  (left  image)  and  the  InSEA  (right  image), 

pipe  was  then  transferred  to  a  mesh  bag  with  3  mm  openings  thus  capturing  only  the 
larger  shells  and  shell  pieces.  The  collected  shell  was  then  dry-sieved  and  analyzed  using 
the  same  procedure  as  was  used  for  the  diver  cores  shells. 

‘Attenuation  Array”  and  InSEA 

To  measure  sound  speed  and  attenuation  in  the  sediment,  two  different  systems  were  used 
over  the  course  of  this  effort.  The  first  is  called  the  “attenuation  array”  and  was 
developed  for  SAX99  and  SAX04  as  part  of  the  larger  Sediment  Transmission 
Measurement  System  (STMSl)  [16].  The  attenuation  array  consists  of  2  sources  and  two 
receivers  mounted  at  the  ends  of  rods  on  a  diver-deployable  frame.  The  transducers  feed 
into  a  junction  box  on  the  frame  and  are  then  cabled  to  the  main  STMSl  control  unit. 
The  diver  pushes  against  the  frame,  driving  the  transducers  into  the  sediment,  and  then 
alerts  the  operator  on  the  ship  that  the  frame  is  in  place.  The  operator  then  triggers  the 
data  collection  and,  once  complete,  asks  the  diver  to  move  the  frame  to  a  new  location. 
During  data  collection,  one  source  transmits  from  40  -  100  kHz  while  the  second  source 
transmits  from  100  -  300  kHz.  The  receivers  are  placed  at  different  distances  from  each 
source  and  the  sound  speed  is  determined  by  the  difference  in  arrival  time  at  each 
receiver.  The  attenuation  is  determined  by  taking  the  ratio  of  the  amplitudes  measured  at 
each  receiver  and  comparing  them  to  the  ratio  measured  when  the  unit  is  in  water  [17]. 

The  Investigation  of  Sediments  by  Acoustics  (InSEA)  system  is  an  improvement  on  the 
design  of  the  attenuation  array  and  was  developed  by  Methodes  Acoustiques  de 
REconnaissance  de  TEnvironnement  in  conjunction  with  Dr.  Laurent  Guillon  at  the 
Erench  Naval  Academy  Research  Institute.  The  system  is  similar  in  design  to  the 
attenuation  array  with  two  sources  and  two  receivers,  each  mounted  in  stakes  that  can  be 
driven  into  the  sediment  by  a  diver.  The  benefit  of  this  system  is  that  all  of  the 
electronics  are  in  a  pressure  case  on  the  frame  and  the  diver  can  initiate  data  collection  by 
pressing  a  button  on  the  case.  The  system  is  very  portable  and  operation  does  not  require 
diver  communications,  as  is  the  case  with  the  attenuation  array.  Eunding  from  ONR- 
Global  supported  the  deployment  of  this  system  during  TREX13. 
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Inversion  of  MBES  data 

The  ultimate  goal  of  any  seabed  classilieation  technique  is  to  determine  the  intrinsic 
physical  properties  of  the  sediment  such  as  density,  sound  speed,  and  porosity.  In  the 
acoustic  detection  of  UXO  on  the  seafloor,  the  signal-to-noise  is  going  to  be  strongly 
affected  by  the  amount  of  scattering  due  to  seafloor  roughness  and  volume  heterogeneity 
within  the  sediment  and  these  become  important  sediment  parameters  in  their  own  right. 
As  opposed  to  the  empirical  approaches  used  in  other  classification  techniques,  in  the 
inversion  developed  here  the  acoustic  data  collected  with  a  MBES  is  used  to  first  obtain 
estimates  of  the  interface  scattering  strength,  volume  scattering,  and  attenuation.  These 
parameters  are  estimated  by  fitting  a  model  for  the  echo  intensity  time  series,  as  in  [7]- 
[10].  This  approach  separates  the  interface  and  volume  scattering,  which  is  possible  due 
to  the  high  directivity  of  the  multibeam  sonar.  The  estimates  of  the  interface  scattering 
strength  can  then  be  used  to  estimate  the  roughness  parameters  and  the  impedance  of  the 
sediment.  Further  estimates  of  the  sediment  parameters  such  as  porosity  and  density  can 
be  obtained  by  using  these  parameters  with  published  empirical  relationships.  A  key 
point  to  this  technique  is  that  all  of  these  parameters  do  not  rely  on  any  additional 
measurements  of  the  sediment  such  as  a  grab  samples  or  in-situ  sound  speed 
measurements.  This  process  is  summarized  in  the  flow  chart  in  Figure  4  and  discussed  in 
detail  below. 


Preprocessing 
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Figure  4:  Inversion  algorithm  flow  chart. 
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The  multibeam  eehosounder  is  typieally  used  to  measure  the  seafloor  bathymetry.  In  this 
application,  the  computer  running  the  sonar  beamforms  the  return  for  each  ping  and  then 
extracts  the  profile  of  the  seafloor  from  the  envelope  of  that  ping.  The  measured  time 
series  data  are  then  either  discarded  or  truncated  about  the  seafloor  return  for  future 
refinements  to  the  bathymetry.  In  order  to  perform  an  inversion  to  extract  the  seafloor 
properties,  the  full  time  series  data  for  each  ping  is  recorded  for  later  processing  in 
MATLAB.  A  full  implementation  of  the  inversion  processing  capable  of  a  wide  area 
survey  would  perform  this  processing  on  the  fly  and,  like  bathymetric  processing,  would 
not  require  the  user  to  retain  the  full,  multi-element  data. 

The  preprocessing  of  the  data  begins,  as  in  the  bathymetric  processing,  with  beamforming 
(using  Hamming  weighting)  with  the  beam  directions  fixed  relative  to  the  sonar  and 
hence  the  ship.  The  time  series  for  each  beam  is  converted  from  complex  to  absolute 
values  (envelope)  and  multiplied  by  calibration  constants.  These  constants  were 
measured  by  both  Teledyne -RESON  and  APL-UW  in  their  respective  acoustic  calibration 
facilities.  With  the  calibration,  the  resulting  time  series  is  the  mean-square  relative 
pressure  for  a  source  giving  unit  RMS  pressure  at  unit  distance. 

Once  the  calibrated  intensity  is  found,  the  pings  are  collected  into  20-ping  sequential 
groups  called  “segments.”  The  20  pings  are  further  subdivided  into  four  “batches”  of 
five  pings  each.  After  aligning  the  leading  edges  of  the  envelope  waveforms,  an  average 
is  taken  to  determine  the  RMS  envelope  of  each  beam  in  each  batch.  Using  these 
averaged  pings,  the  effective  roll  of  the  ship  is  determined  by  looking  across  all  angles 
and  performing  a  linear  fit  to  the  bathymetric  return  as  shown  in  Figure  5.  The  term 


199.995  kHz,  Fit  to  Seafloor  to  Determine  Roll  =  3.1  deg.  Batch  No.  3 
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Figure  5:  Linear  fit  to  the  bathymetric  return  in  the  beamformed  MBES  data.  This 
linear  fit  determines  the  effective  roll  of  the  ship  relative  to  the  seafloor. 
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"effective  roll"  is  used  here  since  roll  is  usually  defined  relative  to  the  true  horizontal.  In 
this  application,  the  goal  of  the  preprocessing  is  to  determine  the  intensity  data  as  a 
function  of  angle  relative  to  the  seafloor.  The  effective  roll  then  determines  the  angular 
offset  of  the  sonar  relative  to  the  cross-track  seafloor  slope.  Once  this  offset  has  been 
applied  to  the  data,  the  batches  are  reoriented  to  the  seafloor  normal  and  the  data  is 
placed  in  angular  bins  that  are  spaced  by  0.5  degrees  from  0  to  50  degrees  off-nadir. 
Port  and  starboard  angles  are  binned  together  and  the  bins  only  refer  to  the  absolute  value 
of  angle  from  nadir.  This  increases  the  data  for  each  scattering  angle  and  simplifies  the 
model-fitting  procedure.  At  this  point,  the  batches  are  averaged  within  each  angular  bin, 
and  the  data  for  each  segment  is  reduced  to  100  intensity  time  series  as  a  function  of 
angle. 

First  Stage 

The  intensity  time  series  generated  in  the  preprocessing  are  then  fit  in  the  first  stage  to 
determine  the  interface  scattering  strength,  volume  scattering  strength,  and  attenuation. 
This  fitting  process  is  explained  in  detail  in  [18],  but  a  brief  summary  is  given  here. 

An  inversion  is  performed  for  each  segment.  The  first  step  in  the  inversion  is  to  fit  a 
simple  time-domain  sonar-equation  model.  The  seafloor  is  assumed  to  be  isotropic  and  is 
represented  by  a  2D  distribution  of  point  interface  scatterers  having  angle-dependent 
interface  scattering  cross  section  and  point  volume  scatterers  with  depth-independent 
volume  scattering  cross  section.  The  attenuation  within  the  sediment  is  assumed  to  be 
depth-independent  as  well.  Time  series  are  generated  separately  for  each  of  the  interface 
scatterers  and  for  all  of  the  volume  scatterers.  The  responses  from  the  interface  and  the 
volume  are  weighted  by  the  initial  guesses  for  the  scattering  strengths  as  well  as  losses 
due  to  effects  such  as  spherical  spreading  and  attenuation  within  the  sediment.  The 
source  and  receiver  directivities  are  accounted  for  in  the  weighting  function  by 
integrating  the  product  of  the  two  beam  patterns  over  the  circle  corresponding  to  the 
grazing  angle  of  interest.  Using  theoretical  beam  patterns,  which  were  found  during 
calibration  to  be  accurate  at  the  beam  angles  of  interest,  a  two-dimensional  look-up  table 
was  generated  assuming  a  fixed  pitch  for  each  field  test  that  allows  rapid  generation  of 
the  time  series.  Figure  6  shows  a  typical  time  series  generated  using  this  approach  where 
the  volume  echo  is  substantially  separated  in  the  time  domain  from  the  interface  echo. 
This  is  due  to  the  delay  in  the  propagation  into  the  seafloor  and  the  effect  of  narrow 
directivity  in  limiting  the  time  duration  of  the  interface  echo.  This  time  separation  is 
critical  to  the  goal  of  quantifying  interface  and  volume  scattering  as  well  as  determining 
attenuation. 

In  fitting,  a  mean-square-error  (MSB)  cost  function  is  used,  and  the  fit  is  made 
simultaneously  to  a  subset  of  the  available  angular  bins  (0  to  30°off-nadir).  This  fit  is 
performed  for  each  20-ping  segment.  This  fitting  process  gives  the  interface  cross- 
section  as  a  function  of  angle,  volume  scattering  cross  section,  and  attenuation  for  each 
segment.  The  MSB  is  quadratic  in  the  two  cross  sections,  so  they  can  be  found  very 
rapidly  by  matrix  inversion  once  a  value  is  chosen  for  attenuation.  The  MSB  is  non¬ 
quadratic  in  attenuation,  so  to  take  advantage  of  the  rapid  inversion  for  cross  sections,  fits 
are  performed  for  10  values  of  loss  parameter,  6,  evenly  spaced  in  the  interval  [0.001 
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time  (s) 

Figure  6:  Example  of  model  time  series  for  interface  (red)  and  volume  (blue) 
scattering,  250  kHz,  8,5°  off-nadir, 

0.01].  The  loss  parameter  giving  the  lowest  MSE  is  chosen.  The  model  does  not  include 
refraction,  and  cannot  determine  volume  scattering  cross  section  directly,  as  it  is  only 
sensitive  to  the  product  of  the  power  transmission  coefficient  and  volume  scattering  cross 
section.  Thus  nominal  sediment  parameters  are  used  to  give  a  rough  value  for  the 
transmission  coefficient.  As  the  transmission  coefficient  only  varies  by  approximately  0.5 
dB  over  the  expected  range  of  sediment  types  (from  mud  to  sand),  the  maximum  error 
incurred  is  0.5  dB  when  the  nominal  parameters  are  badly  chosen.  It  is  possible  to  refine 
this  estimate  after  later  stages  in  the  inversion  are  competed.  Interface  scattering  vs. 
angle  is  represented  by  a  staircase  function  with  step  widths  of  2.5°  out  to  30°.  Data  at 
shallower  angles  are  not  used,  as  the  inversion  becomes  unstable.  This  may  be  due  to  a 
number  of  factors,  such  as  neglect  of  refraction  and  ambiguity  between  interface  and 
volume  scattering  at  shallower  angles. 

At  the  end  of  the  first  stage,  after  fitting,  values  are  obtained  for  each  of  the  following 
parameters: 

1.  Interface  scattering  cross  section  vs.  angle  to  30°  off-nadir 

2.  Volume  scattering  strength  (depth  independent) 

3.  Loss  parameter  {6,  depth  independent) 

Second  Stage 

In  the  first  stage  no  physics  is  employed  beyond  that  of  the  simple  model,  which  is  a 
sonar  equation.  In  the  next  stage,  the  small-slope  (SS)  model  (p.  334  in  [11])  is  fitted  to 
the  interface  scattering  strength  found  in  the  first  stage.  As  SS  requires  sediment  acoustic 
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parameters  as  well  as  roughness  parameters,  nominal  values  are  used  for  density,  sound- 
speed  ratio,  and  loss  parameter.  The  roughness  spectrum  is  taken  to  be  a  pure  power  law 
with  infinite  outer  scale. 


Wik) 


W2 

ky2’ 


(2) 


where  k  is  the  wavenumber,  and  a  global  search  is  performed  for  the  spectral  strength,  W2, 
and  the  spectral  exponent,  72,  to  get  the  best  fit  to  the  angular  shape  of  the  scattering  cross 
section  while  also  determining  the  multiplicative  factor  needed  to  get  a  good  fit.  An 
example  of  this  fit  is  shown  in  Figure  7.  The  multiplicative  factor  is  used  to  correct  the 
nominal  impedance  assumed  at  the  outset.  The  corrected  impedance  is  converted  to 
Richardson’s  Index  of  Impedance,  lOI  =  pvp  where  p  is  the  sediment  mass  density  in 
g/cm^  and  Vp  is  the  sound-speed  ratio. 


Grazing  Angle  (deg) 


Figure  7:  Example  of  fit  of  small-slope  model  to  interface  scattering  strength 
obtained  from  inversion.  This  is  the  sand  section  of  the  TREX13  track;  the 

frequency  is  300  kHz, 


At  the  end  of  the  second  stage,  values  are  obtained  for  each  of  the  following  parameters: 

1 .  Roughness  spectral  strength,  W2 

2.  Roughness  spectral  exponent,  72 

3 .  Index  of  impedance,  lOI 
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Third  Stage 

Using  the  lOI,  the  third  stage  utilizes  empirical  relationships  between  the  physical 
properties  of  the  sediment  and  the  lOI  derived  from  extensive  measurements  made  on 
siliciclastic  and  carbonate  sediments  at  400  kHz.  The  relationships  are  expressed  as 
regressions  in  Table  5.6  of  [6]  and  are  discussed  extensively  there.  Since  these 
regressions  utilize  the  lOI,  which  is  the  product  of  the  sound  speed  and  density,  the  errors 
associated  with  these  regressions  are  small  (r^  =  0.96  -  0.99).  Since  for  these  sediments, 
the  sound  speed  is  strongly  correlated  to  the  porosity,  the  errors  associated  with  the 
porosity  are  also  small  (r^  =  0.99).  The  loss  parameter,  mean  grain  size,  and  sand/gravel 
ratio  have  large  errors  associated  with  them  and  these  are  expected  to  have  the  greatest 
uncertainty  associated  with  them  in  the  final  output. 

The  loss  parameter  obtained  at  this  stage  is  redundant  in  that  the  loss  parameter  was 
already  obtained  in  the  first  stage.  The  application  to  TREX13  and  GulfExl4  indicates 
that  the  loss  parameter  from  the  first  stage  is  smaller  than  the  loss  parameter  from  the 
third,  perhaps  because  the  latter  includes  scattering  effects  and  the  former  is  closer  to  the 
intrinsic  value. 

Using  these  lOI  regressions,  at  the  end  of  the  third  stage,  values  are  obtained  for  each  of 
the  following  parameters: 

1 .  Sound-speed  ratio,  Vp 

2.  Density,  p 

3.  Loss  parameter,  6 

4.  Porosity,  f 

5.  Mean  grain  size 

6.  Sand/Gravel  ratio 

Data  vector 

All  of  the  outputs  of  the  inversion  process  are  summarized  in  the  right  hand  side  of 
Figure  4.  In  all  there  are  12  output  parameters,  only  5  of  which  are  independent  and  thus 
suited  for  classification: 

1.  Volume  scattering  strength 

2.  Loss  parameter  from  stage  1 

3.  Roughness  spectral  strength 

4.  Roughness  spectral  exponent 

5.  Index  of  impedance 

These  outputs  are  calculated  for  each  segment  and  at  each  frequency. 

Additional  outputs 

In  developing  the  inversion  and  understanding  its  performance  in  different  environments, 
several  additional  outputs  have  been  examined  and  could  potentially  be  integrated  into  a 
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final  implementation.  These  outputs  have  to  this  point  been  used  only  qualitatively  or  to 
assess  the  performance  of  the  inversion  algorithm  and  hence  have  not  been  evaluated 
with  the  same  rigor.  They  are  presented  here  as  possible  future  additions  to  this  sediment 
classification  scheme. 

1 .  An  additional  output  that  is  likely  to  prove  useful  is  a  2D  echo-sounder  plot  made 
using  the  nadir  bin.  This  will  be  helpful  in  assessing  layering  and  relative 
strength  of  scattering  by  interfaces. 

2.  Another  output  is  conventional  scattering  cross  section  vs.  angle,  computed  using 
the  sonar  equation  under  the  assumption  that  all  scattering  is  from  the  interface. 
Rather  than  using  a  simple  ct/2  treatment  of  ensonified  area,  where  c  is  the  speed 
of  sound  in  water  and  x  is  the  pulse  width,  the  beam  pattern  azimuthal  integrals 
employed  in  the  more  detailed  model  are  used  to  avoid  problems  at  near-nadir 
angles.  This  type  of  processing  can  be  carried  out  to  shallower  angles  than  the 
inversions  (to  50°  off-nadir  as  of  this  writing).  Penetration  is  so  slight  at  the 
frequencies  of  interest  that  this  scattering  cross  section  agrees  well  with  that 
obtained  by  combining  the  interface  and  volume  contributions  from  the  inversion. 
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Results  and  Discussion 

In  order  to  test  and  evaluate  the  inversion  algorithm,  three  field  tests  were  eondueted;  the 
Gulf  Experiment  in  2011  (GulfExll),  the  Target  and  Reverberation  Experiment  in  2013 
(TREX13),  and  the  St.  Andrew’s  Bay  Experiment  in  2014  (BayExl4).  In  eaeh  ease,  the 
work  presented  here  took  advantage  of  separate  ONR  and  SERDP  funded  experiments  to 
utilize  the  R/V  Sharp  for  the  deployment  of  the  multibeam  sonar  and/or  the  systems  used 
for  the  environmental  eharaeterization  of  the  site.  Also  in  eaeh  ease.  Dr.  Gorm 
Wendelboe,  a  seientist  with  Teledyne -RESON,  partieipated  in  the  experiments  providing 
both  teehnieal  assistanee  and  the  multibeam  sonar  used  to  colleet  the  data  neeessary  to 
evaluate  the  inversion. 

Gulf  Experiment  2011  (GulfExll) 

In  order  to  begin  to  develop  and  test  the  inversion  algorithm,  both  acoustic  and 
environmental  data  were  collected  in  the  Gulf  of  Mexico  in  the  spring  of  201 1  prior  to  the 
start  of  this  project.  These  measurements  were  collected  while  the  R/V  Sharp  was  in  a 
four-point  moor  approximately  two  miles  from  Shell  Island  in  Panama  City  Beach,  PL 
during  an  ONR-sponsored  pilot  study  for  TREX13.  The  experiment  site  had  a  sandy 
seafloor  and  a  water  depth  of  20  m. 

While  the  R/V  Sharp  is  outfitted  with  a  MBES,  Teledyne-RESON  provided  a  calibrated 
Seabat  7125  multibeam  sonar  that  was  mounted  in  the  well  of  the  ship.  The  sonar  was 
modified  for  operation  at  frequencies  of  150,  200,  250,  300,  350,  400,  and  450  kHz  [19]. 
It  employed  a  source  having  a  narrow  beamwidth  in  the  along-track  direction  (1°  at  400 
kHz)  and  a  cross-track  beamwidth  sufficient  to  cover  the  sector  over  which  receiver 
beams  were  formed.  The  near-nadir  receiver  beams  had  a  0.5°  cross-track  width  at  400 
kHz.  Por  frequencies  of  250  kHz  and  lower,  128  beams  were  formed  (from  complex, 
element-level  time  series)  to  cover  a  cross-track  sector  of  ±64°.  At  300  kHz  and  above, 
256  beams  were  formed  to  cover  this  same  sector.  The  sonar  parameters  were  set  to 
suppress  time  varying  gain,  the  ping  rate  was  10  Hz,  and  pulse  lengths  were  160,  120,  96, 
80,  69,  60  and  54  ps  at  the  frequencies  given  above.  Data  were  obtained  in  runs 
comprising  450  to  900  pings,  with  10  to  15  such  runs  at  each  frequency.  One  objective 
was  to  obtain  a  large  number  of  independent  seafloor  echoes  for  subsequent  averaging 
and  comparison  with  models.  Although  adjacent  pings  were  quite  similar  in  envelope 
shape,  ship  movement  inside  the  moorings  during  each  run  was  sufficient  to  provide  the 
desired  statistical  population. 

To  test  the  inversion,  several  measurements  were  made  of  the  environment  at  the 
experiment  site.  At  the  beginning  of  the  experiment,  the  SLS  was  deployed  from  the  R/V 
Sharp  (Pigure  8(a))  in  order  to  measure  the  roughness  of  the  seafloor.  While  the  SLS  is 
mounted  on  the  IMP2,  the  conductivity  probe  was  not  working  at  the  time  of  the  2011 
experiment.  The  motorized  trolley  that  positions  the  conductivity  probe  was  working 
which  allowed  the  SLS  to  collect  a  digital  elevation  map  (DEM)  over  a  3.5  m  long  and  30 
cm  wide  area  on  the  seafloor  with  1  mm  horizontal  resolution.  A  portion  of  the  surface 
measured  during  the  experiment  in  2011  is  shown  in  Pigure  8(b)  with  the  full  DEM 
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Figure  8:  (a)  IMP2/SLS  being  deployed  during  GulfExll,  (b)  Example  of  a  seafloor 
DEM  measured  during  GulfExll  with  tbe  laser  line  intensity  superimposed  to 
distinguish  between  shells  and  sand,  (c)  Complete  DEM  measured  during  a  single 

SLS  deployment  at  the  experiment  site, 

shown  in  Figure  8(c).  While  the  attenuation  array  was  also  deployed  at  this  site,  issues 
with  the  electronics  prevented  this  data  from  being  analyzed. 
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Figure  9:  (a)  Scattering  strengths  at  normal  incidence  determined  from  the 
inversion  (black  lines)  and  the  environmental  characterization  (red  lines).  For  both 
the  red  and  black  lines,  the  interface  scattering  strengths  are  shown  as  the  dot- 
dashed  lines  and  the  volume  scattering  contribution  is  shown  as  the  dashed  lines, 
(b)  Scattering  as  a  function  of  grazing  angle  for  each  frequency. 
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Although  there  was  limited  environmental  data  eollected  at  the  GulfExll  site,  the 
roughness  and  acoustic  data  were  sufficient  to  begin  to  develop  and  evaluate  the 
algorithm.  In  order  to  judge  the  repeatability  of  the  inversion  process  and  to  gain  a  rough 
estimate  of  statistical  error,  data  from  three  MBES  runs  at  each  frequency  were  used  in 
three  separate  inversions.  The  resulting  scattering  cross  sections  were  averaged  and  the 
average  was  converted  to  dB.  Normal-incidence  scattering  strengths  are  shown  in  Figure 
9(a).  The  error  bars  were  computed  by  combining  the  standard  deviation  of  scattering 
cross  section  from  the  three  runs  with  a  systematic  error  of  2  dB,  obtained  by  a  somewhat 
subjective  estimation  of  errors  arising  at  each  stage  in  processing  and  fitting.  There  is  a 
slight  upward  trend  with  frequency,  though  this  trend  does  not  appear  to  be  statistically 
significant.  The  interface  and  volume  contributions  are  similar,  but  scattering  at  the 
interface  appears  to  be  slightly  stronger.  As  noted  above,  the  volume  contribution  is 
comprised  of  contributions  from  two  layers,  but  the  final  results  gave  no  indication  that 
layering  is  important,  that  is,  the  volume  scattering  strengths  of  the  upper  layer  and  the 
remainder  are  comparable. 

We  can  estimate  the  scattering  due  to  surface  roughness  at  the  experiment  site  [20]. 
Figure  10(a)  shows  the  roughness  power  spectrum  determined  from  the  DFM  shown  in 
Figure  8.  Using  the  power  spectral  fit,  the  Kirchhoff  approximation  can  be  used  to  obtain 
the  red  dot-dashed  curve  in  Figure  9(a).  The  difference  between  the  predicted  surface 
scattering  and  the  inversion  results  is  on  the  order  of  10  to  20  dB.  This  discrepancy  was 
initially  thought  to  be  due  to  an  issue  with  the  roughness  measurement.  While  the  SFS 
has  been  extensively  tested  and  undergone  significant  calibration,  there  is  some 
uncertainty  as  to  whether  the  surface  measured  by  the  SFS  was  the  same  surface  seen  by 
the  multibeam.  The  SFS  was  deployed  at  the  beginning  of  the  experiment  and,  after 
several  days,  it  accumulated  a  significant  number  of  fish.  This  fish  density  reached  the 
point  where  it  not  possible  to  collect  roughness  data  due  to  the  occlusion  of  the  camera  by 
the  fish  (Figure  11).  As  a  result,  the  only  useable  roughness  data  was  collected  at  the 
beginning  of  the  experiment.  The  multibeam  data  was  collected  later  in  the  experiment 
and  by  that  time  the  fish  that  were  attracted  to  the  site  most  likely  modified  the  roughness 
spectra.  This  was  thought  to  account  for  the  improvement  in  the  scattering  comparison 
for  angles  away  from  normal  incidence,  as  shown  in  Figure  9(b),  where  the  multibeam  is 
scattering  from  areas  of  the  seafloor  where  the  fish  did  not  rework  the  sediment. 
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(a)  (b) 

Figure  10:  (a)  Roughness  power  spectra  determined  from  the  DEM  shown  in 
Figure  8(c),  The  red  line  is  a  power  law  fit  to  the  roughness  spectra  used  to 
determine  the  two-dimensional  power-law  spectral  fit,  the  parameters  of  which  are 
shown,  (h)  The  shell  size  distribution  measured  at  the  GulfExll  site  (red  and  blue 
circles)  compared  to  the  sand  and  shell  distribution  measured  during  SAX04, 

Subsequent  analysis  of  the  multibeam  data,  however,  indieates  that  the  discrepaney  may 
not  be  due  to  fish-induced  roughness,  but  rather  due  to  issues  related  to  saturation  of  the 
multibeam  data.  Efforts  to  compensate  for  the  saturated  data  may  also  be  responsible  for 
the  areas  of  high  scattering  seen  at  angles  far  from  nadir.  These  saturation  issues  are  now 
well  understood  and  were  avoided  during  TREX13.  While  fish  may  not  have  been 
responsible  for  the  scattering  differences  observed  in  this  experiment,  the  presence  of  fish 
and  their  effect  on  the  seafloor  roughness  remained  an  important  consideration  in 
planning  for  the  2013  experiment. 


While  the  roughness  scattering  data/model  comparisons  may  be  affected  by  the  saturation 
of  the  multibeam  data,  measurements  of  the  shell  size  distribution  at  the  experiment  site, 
shown  in  Figure  10(b)  indicate  that  shell  pieces  make  up  3%  of  the  sediment  volume. 
This  is  significantly  more  than  has  been  found  in  typical  sand  sediments  which  are  on  the 
order  of  1%  [21].  It  is  reasonable  to  assume  that  the  majority  of  scattering  from  the 
sediment  volume  is  due  to  scattering 
from  the  shell  pieces.  By  treating  the 
shell  pieces  as  Rayleigh  or  geometric 
scatterers,  the  volume  scattering 
strength  can  be  estimated  using  the 
shell  size  distribution  shown  in  Figure 
10(b).  The  predicted  scattering 
strength  is  shown  as  the  dashed  red 
curve  in  Figure  9(a).  The  comparison 
between  the  inverted  and  predicted 
values  are  much  better  for  volume 
scattering  than  they  are  for  roughness 
scattering  with  minor  discrepancies  at 

the  lower  frequencies.  These  Figure  11:  IMP2  and  SLS  late  in  GulfExll, 
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discrepancies  may  be  due  to  the 
presence  of  scattering  from 
continuous  volume  heterogeneities 
that  are  not  aeeounted  for  in  the 
model.  The  diserete  scattering 
inversion  model  also  treats  the 
shells  in  the  volume  and  on  the 
surface  equally.  The  shells  on  the 
surface  however  may  have  both  a 
different  size  distribution  and  a 
seattering  response  than  those  in 
the  sediment. 


150 


350 


Frequency  (kHz) 

Figure  12:  Attenuation  at  the  GulfExll  site 
found  from  an  inversion  of  the  multiheam  data. 


As  mentioned  above,  attenuation  is 
obtained  as  part  of  the  fitting 
proeess.  The  results  shown  in 
Figure  12  are  eonsistent  with 
typieal  values  for  sand.  Apart  from 

the  dip  at  200  kHz,  the  data  are  consistent  with  eonstant  attenuation  in  dB/wavelength. 
This  is  equivalent  to  a  linear  rise  in  attenuation  with  frequency  that  is  consistent  with  the 
frequeney  dependenee  measured  in  the  laboratory  and  in  other  sand  sediments.  Due  to 
problems  with  the  attenuation  array,  attenuation  was  not  measured  at  the  experiment  site 
and  cannot  be  compared  with  the  inversion  results. 


Target  and  Reverberation  Experiment  2013  (TREX13) 


The  second  field  test  took  place  in  the  spring  of  2013  as  part  of  TREX13.  The  main 
experiment  was  composed  of  two  parts:  an  ONR  sponsored  reverberation  experiment  and 
a  synthetic  aperture  sonar  (SAS)  target  deteetion  and  elassification  experiment  sponsored 
jointly  by  ONR  and  SERDP.  The  loeation  and  layout  of  the  experiment  is  shown  in 
Figure  13.  A  majority  of  the  assets  shown  in  the  figure  were  deployed  for  the 
reverberation  portion  of  the  experiment  while  the  SAS  work  took  plaee  at  the  moored 
location  of  the  RA^  Sharp  at  the  interseetion  of  the  elutter  and  main  tracks.  The  ship  was 
in  a  four-point  moor  for  a  majority  of  the  experiment,  with  the  exception  of  the  last  two 
weeks  when  it  was  used  to  colleet  environmental  and  bathymetric  data  along  the  main 
traek.  It  was  during  this  time  that  we  were  able  to  utilize  the  ship  to  colleet 
environmental  and  acoustic  data  along  a  350  m  line  to  the  southeast  of  the  mooring 
loeation.  This  line  was  chosen  using  the  bathymetric  and  backscatter  data  collected  by 
Dr.  Christian  de  Moustier  and  Dr.  Barbara  Kraft  during  the  survey  cruise  that  took  plaee 
immediately  before  the  RA^  Sharp  was  moored.  As  shown  in  Figure  14,  the  line  starts  at 
near  the  top  of  a  ridge  and  traverses  a  swale  as  it  moves  towards  the  top  of  the  next  ridge. 
From  the  baekseatter  data,  the  ridge  tops  were  areas  of  high  baekscatter  while  the  ridge 
swales  were  areas  of  low  backscatter.  At  the  time  it  was  believed  that  this  indicated  sand 
along  the  ridges  and  mud  in  swales  and  this  location  allowed  us  to  eolleet  data  in  both 
areas.  This  site  was  also  far  enough  away  from  the  mooring  site  so  as  not  to  be  affeeted 
by  the  fish  that  had  aeeumulated  around  the  ship  and  the  deployed  instruments. 
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Figure  13:  Bathymetry  and  distribution  of  assets  for  main  reverberation  component 
of  TREX13,  The  location  of  the  MBES  Study  Site  to  support  the  multibeam 
inversion  work  indicated  to  the  southeast  of  the  R/V  Sharp, 


(a)  (b) 

Figure  14:  (a)  Backscatter  measured  using  the  multibeam  sonar  during  the  survey 
cruise  at  the  start  of  TREX13,  The  multibeam  measurement  track  is  shown  as  the 
red  hue,  (b)  The  scattering  strength  (upper  panel)  and  bathymetry  (lower  panel) 

along  the  inversion  line. 

The  collection  of  MBES  and  environmental  data  took  place  from  May  25  to  June  4,  2013. 
The  IMP2  and  SLS  was  deployed  along  the  track  on  May  25,  multibeam  data  was 
collected  along  the  track  on  May  27  and  28,  and  divers  were  deployed  along  the  track  on 
June  2  and  June  4  to  collect  cores,  samples,  and  deploy  the  InSEA.  There  were  no 
weather  events  and  no  significant  fish  accumulation  during  this  time. 

MBES  Data 

Eor  TREX13,  Dr.  Wendelboe  brought  two  Seabat  7125s  from  Teledyne-RESON;  the 
original  system  used  during  GulfExll  and  a  newer  prototype  that  had  several 
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modifications  that  helped  to  avoid  the 
saturation  issues  encountered  during 
GultExll.  Data  was  collected  using  both 
systems,  but  the  latter  system  produced 
higher  quality  data  that  was  eventually  the 
focus  of  the  inversion  testing.  Due  to  time 
limitations  and  to  facilitate  the  deployments 
of  both  sonars,  the  7125s  were  deployed 
using  a  pole  mount  on  the  side  of  the  ship 
(Figure  15).  While  this  enabled  us  to 
change  out  the  sonars,  this  deployment 
configuration  is  unstable  and  can  experience 
significant  vibrations  while  the  ship  is 
underway.  In  an  attempt  to  mitigate  this 
issue,  the  ship  used  dynamic  positioning  to 
crab  along  the  350  m  track  at  slow  speeds 
(~1  kn).  Also,  since  the  sonar  wasn’t 
exactly  aligned  with  the  ship  as  would  be 
the  case  in  a  well  mount,  there  were 
unknown  offsets  in  the  pitch  and  roll  of  the 
sonar.  The  preprocessing  for  the  inversion 
extracts  the  roll  information  so  this  isn’t  a 
problem,  but  the  pitch  does  affect  the 
inversion  performance  and  was  estimated 
through  trial  and  error  to  be  approximately  3° 


Like  the  GulfExll  sonar,  the  new  sonar  was  also  capable  of  operating  at  multiple 
frequencies  and  during  this  experiment  both  systems  were  operated  at  frequencies  of  180, 
200,  225,  250,  275,  300,  325,  350,  375,  400,  and  420  klTz.  For  each  frequency,  the  RW 
Sharp  started  at  one  end  of  the  track  and  crabbed  slowly  to  the  other  end  where  the 
frequency  was  changed  for  the  return  trip.  This  technique  is  not  amenable  for  a  wide- 
area  survey  such  as  the  bathymetric  survey  shown  in  Figure  13  but  worked  well  over  the 
small  track  to  collect  data  to  test  the  inversion  algorithm. 

Environmental  Data 


To  evaluate  the  accuracy  of  the  MBES  inversion,  environmental  measurements  were 
collected  every  50  meters  along  the  track  for  a  total  of  eight  locations  as  shown  in  Figure 
16.  Ideally  every  system  and  sampling  technique  would  be  used  at  each  location,  but 
due  to  logistical  constraints  and  the  limited  time  allotted  for  these  measurements  only  a 
few  sites  were  fully  characterized.  Table  1  lists  each  of  the  environmental  measurements 
and  at  which  locations  they  were  made. 


Figure  15:  Pole  mount  used  deploy  the 
Seabat  7125  sonars  from  the  RA^ 
Sharp  during  TREX13, 
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Figure  16:  Locations  of  the  environmental  measurements  along  the  MBES  track. 


System/Technique 

Environmental 

Measurement 

Sites 

Conductivity  Probe 

Volume  Heterogeneity 

1,5,8 

Seafloor  Easer  Scanner 

Bottom  Roughness 

1,2,3,4,6,7,8 

Diver  Cores 

Bulk  Porosity 

Shell  Size  Distribution 

1, 4,5,8 

Sediment  Samples 

Shell  Size  Distribution 

1, 4,5,8 

In-Sediment  Acoustic  Probe 

Sound  Speed  and  Attenuation 

1, 4,5,8 

Table  1:  Environmental  measurements  and  locations  at  along  the  MBES  track. 


Seafloor  roughness  was  measured  at  all  sites  along  the  track  except  for  site  #5.  In  the 
backscatter  data  shown  in  Figure  17,  site  #5  lies  on  the  edge  of  the  line  of  extremely  low 
backscatter.  This  region  is  mud  that  has  collected  at  the  base  of  the  sand  ridge  and  the 
IMP2’s  pads  sunk  10-15  cm  into  this  mud.  This  caused  the  laser  line  to  be  projected 
outside  the  field  of  view  of  the  camera  and  data  could  not  be  collected.  The  ID 
roughness  power  spectra  measured  at  the  other  locations  along  the  measurement  track  are 
shown  in  Figure  17.  At  the  high  wavenumbers,  the  spectra  measured  on  the  ridge  are  5- 
10  dB  higher  than  the  spectra  measured  in  the  swale.  Examination  of  the  DEMs 
corresponding  to  these  spectra  indicates  this  higher  roughness  is  driven  by  the  shells 
which  cover  the  seafloor  on  the  ridge.  These  shells  are  absent  in  the  swales  and  hence  the 
interface  is  not  as  rough.  Site  #4  has  significantly  higher  roughness  due  to  the  increase  in 
the  number  of  shells  as  the  mud  area  is  approached  from  the  west. 

In  order  to  model  the  scattering  from  the  seafloor  due  to  roughness,  the  ID  roughness 
power  spectrum  at  each  site  is  fit  by  the  sum  of  two  spectra, 

=  M4(/c)  +  W^ndk),  (3) 

where  n  denotes  the  site  number.  For  wavenumbers  less  than  30  m'\  the  spectra  can  be 
approximately  described  by  the  mean  of  all  of  the  spectra  along  the  track.  Fitting  a  von 
Karman  spectrum  to  this  mean  spectrum  yields 
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Figure  17:  ID  roughness  power  spectra  measured  along  the  MBES  measurement 

track. 
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where  Wa  =  0.003  Ja  =  4.23,  and  £«  =  0.1  m.  For  A:  >  30  m'',  the  speetrum  at  eaeh 
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where  Whn  and  ybn  are  determine  from  a  power  law  fit  at  each  site  to  the  portion  of  the 
spectrum  in  the  range  100  m"'  <  k  <  1300  m'\  The  values  from  these  fits  are  shown  in 
Figure  18(a)  and  (b).  The  outer  scale  is  chosen  to  be  the  same  at  each  site,  Lb  =  0.01  m. 
A  comparison  of  this  combined  fit  to  the  spectra  measured  at  site  #1  and  #8  is  shown  in 
Figure  18(c).  The  2D  roughness  power  spectrum  is  required  to  calculate  seafloor 
scattering  using  the  small-slope  approximation.  Since  the  DEMs  measured  along  the 
MBES  measurement  indicate  that  for  the  high  wavenumbers  the  seafloor  is  isotropic,  the 
spectral  parameters  for  the  2D  power  spectra  can  be  recovered  from  the  ID  spectral 
parameters  using  Eq.  (D.18)  and  (D.19)  in  [1 1]. 

The  deployment  of  the  SES  along  the  track  took  place  over  the  course  of  a  day  and  at 
several  of  the  sites  during  its  deployment,  the  conductivity  probe  was  used  to  measure  the 
porosity  fluctuations  in  the  sediment.  While  the  variance  of  the  fluctuations  was  found  to 
be  higher  at  site  #1  than  at  site  #2,  in  both  cases  scattering  from  these  heterogeneities  was 
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Figure  18:  (a)  Spectral  exponent  and  (b)  spectral  strength  determined  from  a  power 
law  fit  to  the  roughness  power  spectra  measured  along  the  MBES  measurement 
track,  (c)  Roughness  power  spectra  and  the  spectral  fits  for  sites  #1  and  #8, 


Site 

Density  (g/cm^) 

Porosity 

B  (10'^  1/m^) 

P 

1 

1.89 

0.35 

1.539 

3.75 

4 

1.78 

0.33 

1.054 

3.32 

8 

1.96 

0.40 

0.270 

4.07 

Table  2:  Density  and  porosity  determined  from  diver  cores  collected  along  the 

multibeam  track. 


predicted  to  be  significantly  lower  than  that  due  to  scattering  from  the  shell  and  shell 
hash  within  the  sediment.  As  a  consequence,  this  data  has  not  been  used  in  the 
inversion/model  comparisons  discussed  below.  At  site  #5,  although  the  IMP2  footings 
sunk  into  the  mud,  it  was  still  possible  to  deploy  the  conductivity  probe  and  a  12  cm  deep 
and  20  cm  long  set  of  data  was  collected  in  the  mud.  Inclusions  of  sand  were  observed  in 
this  data  but  there  is  too  little  data  to  collect  meaningful  statistics  at  this  site. 


Diver  cores  were  collected  at  sites  #1,  #4,  #5,  and  #8.  Each  of  these  cores  was  analyzed 
for  sediment  density,  porosity,  and  shell  distribution.  The  density  and  porosity 
determined  from  the  cores  are  shown  in  Table  2.  Site  #5  is  not  shown  since  the  mud  area 
has  not  been  included  in  the  data/model  comparisons.  The  shell  weight  distribution  and 
the  shell  number  density  were  determined  from  both  the  diver  cores  and  larger 
cofferdams.  The  combined  results  of  these  measurements  are  shown  in  Figure  19.  The 

-5 

density  of  the  shell  material  was  not  measured,  but  a  typical  value  of  2.7  g/cm  is 
assumed.  Both  the  shell  weights  and  numbers  are  consistent  with  the  SLS  measurements 
that  indicate  that  the  amount  of  shell  increases  as  the  mud  area  is  approached  from  the 
west  and  then  decreases  abruptly  at  the  mud  area.  The  number  density  is  the  primary 
input  to  the  scattering  model  used  in  the  data/model  comparisons.  A  power  law  has  been 
fit  to  this  data,  as  shown  in  the  lower  panel  of  Figure  19,  which  can  be  expressed  as 


(6) 
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Figure  19:  Shell  weight  distribution  (Upper  panel)  and  shell  numher  density 
(lower  panel)  for  site  #1  (blue),  site  #4  (red),  and  site  #8  (black). 


Figure  20:  Sound  speed  (upper  panel)  and  loss  parameter  or  wavenumber  ratio 
(lower  panel)  measured  using  the  InSEA  along  the  multibeam  track. 


where  a  is  the  shell  radius  in  mm  and  ao  =  I  mm.  The  values  for  the  number  density 
amplitude,  B,  and  slope,  P,  are  given  in  Table  2. 


In  addition  to  the  diver  eores,  sediment  sound  speed  and  attenuation  were  also  colleeted 
at  sites  #1,  #4,  and  #8.  The  results  of  these  measurements  are  shown  in  Figure  20,  where 
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the  sound  speed  and  loss  parameter  (wavenumber  ratio)  are  shown.  The  loss  parameter  is 
expressed  as 


K 


(7) 


where  k'p  is  the  real  part  of  the  sediment  wavenumber  and  kp  is  the  imaginary  part  of  the 
sediment  wavenumber.  The  loss  parameter  measured  using  the  InSEA  eaptures  both  the 
intrinsic  attenuation  that  would  be  present  in  a  sand  sediment  in  the  absence  of  shells  and 
the  attenuation  due  to  scattering  from  the  shells  in  the  sediment. 


Inversion  Results  And  Comparisons  With  Measured  Predictions 


To  present  the  output  of  the  inversion  in  a  way  that  is  consistent  with  the  processing  and 
useful  to  the  user,  the  groups  of  segments  are  associated  with  patches  of  the  seafloor 
along  the  data  collection  track.  The  inversion  uses  the  beams  from  nadir  out  to  30 
degrees  on  either  side.  These  beams  are  used  to  produce  the  data  vector  discussed  in  the 
previous  section  and  there  is  no  direct  angular  dependence  to  the  inverted  parameters. 
The  data  associated  with  a  given  segment  therefore  corresponds  to  the  portion  of  the 
seafloor  of  the  seafloor  interrogated  by  the  60-degree  group  of  beams  that  for  the 
TREX13  water  depth  corresponds  to  a  22  m  swath  width.  In  order  to  present  the 
inversion  in  a  meaningful  way  and  to  evaluate  the  uncertainties  in  the  data,  we  have 
chosen  to  group  the  segments  into  square  patches  with  the  length  of  each  side 
corresponding  to  the  length  of  the  swath  width.  Data  from  all  of  the  segments  within  this 
patch  are  averaged  and  the  confidence  limits  of  the  data  are  determined.  An  example  of  a 
map  of  the  density  values  produced  by  the  inversion  along  the  TREX13  track  is  shown  in 
Figure  21.  This  type  of  display  (with  multiple  adjacent  tracks)  may  be  a  final  data 
product  and  visualization  for  the  user  in  a  future  implementation  of  the  inversion. 


From  the  MBES  inversion,  there  are  five  independent  parameters  which  ideally  could  be 
compared  to  ground  truth  data  in  order  to  evaluate  the  performance  of  the  inversion 
algorithm.  Unfortunately,  these  outputs  either  do  not  correspond  directly  to  a  measurable 
quantity,  as  in  the  case  of  the  lOI,  or  can  be  difficult  to  compare  on  a  one-to-one  basis,  as 
is  the  case  with  the  roughness  spectral  strength  and  exponent.  For  the  roughness  spectra, 
different  values  of  the  spectral  strength  and  exponent  can  produce  similar  spectra  over  a 
limited  range  of  wavenumbers.  It  is  better  in  this  case  to  compare  the  resulting  spectra 
directly  to  determine  how  accurately  the  roughness  spectrum  has  been  determined  and 
over  what  range  the  results  are  valid.  We  will  therefore  focus  on  7  parameters  for  the 
inversion  evaluation; 
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Eastings  (m) 


Figure  21:  Map  of  the  density  ratio  along  the  TREX13  track  determined  from  the 
MBES  inversion.  The  background  is  the  scattering  strength  measured  hy  de 
Moustier  and  Kraft  during  the  experiment. 


1.  Volume  Scattering  Strength 

2.  Loss  Parameter  from  Stage  1 

3.  Roughness  Power  Spectrum 

4.  Density  Ratio 

5.  Sound  Speed  Ratio 

6.  Loss  Parameter  from  lOI 
1.  Porosity 

Comparisons  of  these  parameters  to  the  ground  truth  measurements  are  shown  in  Figure 
22  -  Figure  26. 

A  majority  of  the  parameters  are  essentially  frequency  independent  with  the  exception  of 
the  volume  scattering  strength,  which  is  shown  as  a  function  of  track  position,  and 
frequency  in  Figure  22.  At  each  frequency  there  are  a  number  of  track  positions  where 
there  is  no  data.  To  collect  the  MBES  data,  recall  that  the  R/V  Sharp  crabbed  back  and 
forth  along  the  track,  collecting  data  at  a  single  frequency  for  each  run  down  the  track. 
Moving  slowly  along  the  track  made  the  ship  susceptible  to  significant  roll  by  the 
incoming  swell  and  this,  coupled  with  the  instability  of  the  mounting,  corrupted  some  of 
the  segments  along  the  track.  A  rejection  criterion  was  added  to  the  inversion  algorithm 
to  remove  these  extreme  events  from  the  data.  To  invert  for  the  frequency- independent 
parameters,  the  algorithm  uses  the  data  at  all  frequencies  and  as  a  result,  the  dropouts 
tend  to  be  absent  in  these  data  products  but  they  do  increase  the  uncertainty  in  the 
estimate. 
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Figure  22:  Volume  scattering  strength  along  the  TREX13  track  (blue)  compared  to 
model  predictions  made  using  the  measured  sediment  properties  and  shell 

distributions  (red). 
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Figure  23:  Volume  scattering  strength  measured  at  0  m,  150  m,  and  350  m  along  the 
TREX13  track  (blue)  compared  to  ground  truth  predictions  (red). 
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The  volume  scattering  data  is  compared  to  predictions  made  using  a  discrete  scattering 
theory  developed  by  Dr.  Anatoliy  Ivakin  [22],  which  uses  as  input  the  measured  shell 
number  density  and  sediment  parameters  (Figure  22  and  Figure  23).  The  volume 
scattering  model  uses  only  the  number  density  values  measured  along  the  track,  which,  as 
can  seen  in  Figure  19,  extend  from  3  mm  to  9.5  mm.  In  reality,  there  are  shells  smaller 
than  3  mm  but  we’ve  chosen  not  to  extrapolate  to  these  smaller  sizes.  The  flat  frequency 
response  seen  in  Figure  23  is  a  consequence  of  this  truncation  of  the  number  density  at  3 
mm.  Extrapolating  to  lower  sizes  would  produce  a  curve  that  increases  with  frequency 
and  would  be  roughly  4  dB  higher  at  200  kHz  and  9  dB  higher  at  400  kHz.  These 
discrepancies  between  the  volume  scattering  strength  and  the  inversion  for  the  portion  of 
the  track  less  than  200  m  may  be  due  to  the  failure  of  the  volume  scattering  model  to 
account  for  multiple  scattering.  There  is  significant  shell  content  at  this  site  and  if 
multiple  scattering  were  included  in  the  model  we  would  expect  the  scattering  strength  to 
decrease. 

The  inversion  calculates  the  loss  parameter  at  two  different  stages  in  the  processing.  The 
first  takes  place  in  the  fitting  procedure  during  the  first  stage  while  the  second  is  derived 
from  the  lOI  that  is  found  in  the  second  stage.  Both  of  these  results  are  shown  as  a 
function  of  distance  along  the  track  in  Figure  24.  The  loss  parameter  from  the  lOI  is 
roughly  two  to  three  times  as  large  as  that  found  in  the  first  stage.  The  value  from  the  lOI 
at  the  beginning  and  end  of  the  track  also  compare  well  to  the  values  measured  using  the 
InSEA.  This  is  consistent  with  the  hypothesis  that  the  loss  parameter  in  the  first  stage 
includes  only  the  intrinsic  attenuation  and  scattering  from  the  shell  pieces  or  other 
inhomogeneities  is  captured  as  volume  scattering.  Eor  the  InSEA,  any  energy  scattered 
from  the  coherent  field  is  lost  in  the  measurement  and  manifests  as  additional  attenuation. 
Since  the  lOI  are  derived  from  similar  types  of  measurements,  it  is  logical  that  the  lOI 
loss  parameter  should  also  include  scattering  loss.  This  picture  is  not  consistent  however 
with  the  increase  in  attenuation  at  the  end  of  the  track  where  the  shell  content  is 
significantly  lower. 

The  roughness  spectra  produced  using  the  spectral  strength  and  exponent  produced  by  the 
inversion  are  shown  in  Eigure  25  as  a  function  of  range  for  k  =  100,  215,  464,  and  1000 
rad/m.  Typically  the  roughness  spectrum  is  shown  as  a  function  wavenumber  as  in 
Eigure  18(c),  but  this  type  of  plot  makes  it  difficult  to  compare  the  inverted  spectrum  to 
the  measured  spectrum.  Over  this  band  of  wavenumbers,  the  inverted  spectra  quite 
clearly  capture  the  transition  from  shelly  sand  to  the  mud/sand  side  and  the  change  in 
spectral  level  is  comparable  to  the  measured  change.  The  absolute  levels  compare  quite 
well  at  the  highest  wavenumbers  while  the  inversion  overestimates  the  spectra  as  the 
wavenumber  decreases.  At  k  =  10  rad/m  (not  shown  in  Eigure  25),  this  overestimation  is 
on  the  order  of  10  dB.  The  cause  of  this  divergence  is  not  clear.  In  applying  these  results 
to  a  sonar  system  such  as  an  SAS,  which  operates  at  shallow  grazing  angles,  a  bragg 
wavenumber  of  100  rad/m  corresponds  to  a  sonar  frequency  of  12  kHz.  At  this  frequency 
we  would  expect  a  roughly  5  dB  error  while  for  the  lower  two  panels  of  Eigure  25,  we 
expect  the  inversion  results  to  apply  quite  well  to  a  sonar  operating  down  to  a  least  50 
kHz. 
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Figure  24:  Loss  Parameter  determined  from  stage  1  of  the  inversion,  from  the 
inversion  lOI  result,  and  the  InSEA, 


0  50  100  150  200  250  300  350 

Range (m) 


Range (m)  Range (m) 


Figure  25:  Roughness  Spectra  along  the  TREX13  track  compared  to  the  spectra 
measured  using  the  SLS,  Each  panel  corresponds  to  a  different  wavenumber. 
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Several  of  the  values  determined  using  the  lOI  from  the  inversion  are  shown  in  Figure  26. 
The  lOI  outputs  compare  reasonably  well  in  the  shelly  sand  area,  slightly  overestimating 
the  values,  but  diverge  significantly  past  the  200  m  mark.  The  goal  of  this  field  test  was 
to  examine  a  sand  sediment  and,  as  a  result,  two  of  the  three  ground  truth  measurements 
were  focused  on  the  shelly  sand  side.  The  inversion  appears  to  be  tending  towards  the 
values  measured  at  350  m  but  there  is  not  enough  information  in  the  200  -  350  m  range  to 
test  this  aspect  of  the  inversion  here. 


Range  (m) 


Figure  26:  Sound  Speed  Ratio  (upper),  Porosity  (middle),  and  Density  Ratio  (lower) 
determined  using  the  inverted  lOI  (blue  points).  The  measured  values  are  also 

shown  (red  points). 
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Despite  the  instability  of  the  MBES  mounting  during  this  field  test,  the  inversion 
performed  admirably  in  a  trying  environment.  The  TREX13  site  is  eertainly  a  diffieult 
one,  much  more  so  than  a  majority  of  sites  visited  in  past  high-frequency  acoustic 
experiments.  Multiple  scattering  seems  to  be  prevalent,  and  the  inversion  model  handles 
this  problem  gracefully.  It  fits  a  single-scattering  model  to  multiple-scattering  data. 
Codes  such  as  PC  SWAT  also  uses  single-scattering  approximations  and,  if  it  uses  the 
inversion  results  as  inputs,  will  produce  the  correct  reverberation  levels  at  our  survey 
frequencies,  if  not  substantially  lower  as  indicated  by  the  spectral  inversions. 

St,  Andrew’s  Bay  Experiment  2014  (BayExl4) 

The  third  field  test  took  place  in  the 
spring  of  2014  as  part  of  BayExl4. 

The  main  experiment  was  again 
focused  on  SAS  target  detection  and 
classification  and  was  jointly 
sponsored  by  ONR  and  SERDP.  The 
location  and  layout  of  the  experiment 
is  shown  in  Eigure  27.  The  RW 
Sharp  was  again  placed  in  a  four- 
point  moor  for  the  entire  experiment 
and  the  environmental  measurements 
supporting  the  MBES  inversion 
evaluation  were  conducted  while  the 
ship  was  moored.  The  location  of 
the  RW  Sharp  was  dictated  by  the 
needs  of  the  main  experiment  and  the 
multibeam  data  track  was  chosen  to 
pass  along  the  south  side  of  the  R/V 
Sharp  location.  The  track  was  200  m 
in  length  with  the  ship  located  at 
roughly  the  100  m  distance. 

The  collection  of  MBES  data  took  place  prior  to  the  beginning  of  the  experiment  on 
March  30  and  April  1,  2014  in  order  to  collect  data  from  the  site  prior  to  the  deployment 
of  equipment  and  the  accumulation  of  fish.  The  ship  entered  the  4-point  moor  on  May  8 
and  remained  there  until  May  27. 

MBES  Data 

For  BayExl4,  Dr.  Wendelboe  provided  three  Seabat  7125s;  the  original  system  used 
during  GulfExll,  the  system  used  during  TREX13,  and  a  new  system  that  included 
design  improvements  over  the  TREX13  system.  Data  was  collected  along  the  MBES 
track  with  each  system.  The  TREX13  sonar  was  a  prototype  in  which  Teledyne-RESON 
had  been  experimenting  with  new  manufacturing  techniques.  As  a  result,  the  sonar 
performance  had  degraded  between  TREX13  and  BayExl4  and  we  used  the  new  system 


Figure  27:  Configuration  of  the  RW  Sharp 
and  mooring  lines  relative  to  the  MBES  track 
The  inset  shows  the  mooring  location  relative 
to  the  shoreline  along  the  hay. 
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Figure  28:  (Left  panel)  Formation  factor  measured  using  the  IMP2  at  the  BayExl4 
site,  (Right  panel)  The  mean  spectra  for  the  mud/water  and  sand/mud  interfaces 
extracted  from  the  formation  factor  data. 


for  the  inversion  evaluation  diseussed  below.  Instead  of  using  the  RA^  Sharp  to  colleet 
the  MBES  data,  a  dedieated  survey  boat  was  used  to  avoid  the  logistical  issues  that  come 
with  the  use  of  the  ship  and  to  facilitate  the  deployment  of  the  multiple  sonars.  The 
workboat  was  hired  from  Seaside  Engineering  and  Surveying,  EEC  and  was  a  26-foot 
aluminum  boat  with  a  moon  well  for  deploying  the  sonars.  The  bay  was  very  calm 
during  the  data  collection  and,  with  the  sonars  mounted  in  the  moon  well,  the  data  was  of 
much  higher  quality  than  TREX13. 

Environmental  Data 


The  sediment  in  St.  Andrew’s  Bay  at  the  BayExl4  site  consisted  of  a  mud  layer  over  a 
sand/mud  sediment.  The  mean  thickness  of  the  layer  was  determined  from  conductivity 
probe  measurements  to  be  13.04  ±  1.25  cm  thick.  An  example  of  a  conductivity  probe 
measurement  at  the  site  is  shown  in  Eigure  28.  To  measure  the  roughness  of  both  layers, 
the  formation  factor  data  from  the  conductivity  probe  was  processed  to  extract  the  ID 
profiles  of  both  interfaces.  The  mean  spectra  are  shown  in  Figure  28  and  the  ID 
roughness  spectral  strength  and  exponents  are  given  in  Table  3.  The  roughness  spectrum 
for  the  mud/water  interface  agreed  with  the  SLS  measurements  and  the  laser  data  showed 
that  the  power  law  continued  to  higher  wavenumbers  (<1000  rad/m).  The  sand/mud 
interface  is  very  rough  and  is  5-10  dB  rougher  than  the  sediment  at  TREX13.  This 
interface  is  also  exceptional  in  that  it  is  very  difficult  to  determine  exactly  where  the 
interface  is  located.  The  surface  shown  in  Figure  28  is  one  an  attempt  to  define  the 
interface  and  alternative  interpretation  would  be  to  define  a  mean  surface  and  interpret 
the  variations  not  as  roughness  but  instead  as  volume  heterogeneity.  As  will  be  discussed 
below,  this  may  explain  why  the  inversion  overestimates  the  volume  scattering  expected 
for  regions  below  the  troughs  of  the  mud-sand  interface. 


ID  Spectral  Strength  (w;) 

ID  Spectral  Exponent  (y) 

Mud/Water  Interface 

0.00014  m'-’' 

2.49 

Sand/Mud  Interface 

0.00075  m'-’’ 

2.06 

Table  3:  Spectral  strengths  and  exponents  for  the  BayExl4  sediment  interfaces. 
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Figure  29:  Density  fluctuation  power  spectra  within  the  underlying  sand  sediment. 


Frequency  (kHz)  Frequency  (kHz) 

Figure  30:  Sound  speed  ratio  and  Attenuation  measured  in  the  mud  layer  (left 
panels)  and  the  underlying  sand  (right  panel). 
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Diver  cores  were  again  used  to  measure  the  sediment  properties,  both  in  the  mud  layer 
and  the  underlying  sand.  There  was  no  significant  shell  content  in  the  underlying  sand 
and  the  shell  size  distribution  was  not  measured.  The  mud  had  a  density  of  y9  =  1 .26  ± 
0.06  g/cm^  and  a  porosity  of  =  81%  ±  1%.  The  underlying  sand  had  a  density  of  p  = 
2.10  ±  0.02  g/cm^  and  a  porosity  ofy5  =  40%  ±  1%. 

Using  the  diver  core  data,  the  formation  factor  data  from  the  conductivity  probe  could  be 
inverted  to  get  the  density  distribution  in  the  mud  layer  and  the  underlying  sand.  The 
data  was  used  to  calculate  the  density  fluctuation  power  spectrum  shown  in  Figure  29. 
There  were  three  different  IMP2  data  sets  used  to  calculate  the  density  fluctuation  power 
spectra  and  the  power  law  fits  to  each  run  yield  the  spectral  parameters  in  Table  4.  These 
measurements  do  not  account  for  heterogeneities  near  the  sand/mud  interface  which  may 
be  significant. 


Spectral  Strength 

Spectral  Exponent 

Run  299 

0.0029  m*-’’ 

1.75 

Run  316 

0.0031  m*-’’ 

1.66 

Run  309 

0.0004  m*-’' 

1.50 

Table  4:  Spectral  strengths  and  exponents  for  the  density  fluctuation  power  spectra 

within  the  sand  sediment. 


To  measure  the  sound  speed  and  attenuation  in  both  the  mud  and  sand  layers,  the 
attenuation  array  was  deployed  from  the  R/V  Sharp.  In  order  to  measure  the  sound  speed 
and  attenuation  in  the  mud  layer,  the  array  had  to  be  reconfigured  bring  the  receivers  very 
close  to  the  source.  This  greatly  increased  the  uncertainty  in  the  measurements  as  can  be 
seen  in  Figure  30.  In  the  mud,  the  mean  sound  speed  ratio  across  the  band  was  0.984  ± 
0.071  and  the  mean  attenuation  was  36.3  ±5.3  dB/m.  For  the  sand  measurements,  one  of 
the  receivers  was  corrupted  by  cross  talk  and  the  processing  had  to  be  done  using  a  single 
receiver.  This  introduced  significant  uncertainty  to  the  sound  speed  ratio  as  can  be  seen 
in  Figure  30  where  the  mean  sound  speed  ratio  is  1.075  ±  0.042  and  a  potentially  large 
systematic  bias  to  attenuation  measurement,  which  has  a  mean  of  46.1  ±3.8  dB/m.  This 
value  is  potentially  an  order  of  magnitude  too  low. 

Inversion  Results  and  Comparisons  With  Measured  Parameters 

For  the  BayExl4  inversion,  we  will  again  focus  on  same  the  7  parameters  examined  for 
the  TREX13  inversion.  Comparisons  of  these  parameters  to  the  ground  truth 
measurements  collected  at  the  R/V  Sharp  position  are  shown  in  Figure  31-  Figure  36. 
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Figure  31:  The  inversion  volume  scattering  strength  along  the  BayExl4  track 
(blue)  compared  to  scattering  strength  determined  from  the  measured  density 
fluctuations  (red)  in  the  sand  sediment. 


Figure  32:  Inverted  volume  scattering  strength  (blue)  at  the  100  m  distance  along 
the  BayExl4  track  compared  to  the  scattering  predictions  (red). 
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As  opposed  to  the  TREX13  inversions,  the  inversions  for  BayExl4  are  mueh  more  stable; 
there  are  very  few  dropouts  in  the  volume  scattering  strength  along  the  track  (Figure  31). 
The  scattering  strengths  are  also  quite  stable  with  only  small  variations  (2-5  dB)  along  the 
track.  The  inversion  seems  to  indicate  that  the  ends  of  the  track  have  a  slightly  higher 
volume  scattering  and  roughness  scattering  response  as  can  be  seen  in  both  Figure  3 1  and 
Figure  35.  The  volume  scattering  at  the  R/V  Sharp  site  was  calculated  using  small- 
perturbation  theory  (pp.  383-385  of  [11])  using  the  sediment  parameters  and  the 
fluctuation  spectra  in  Figure  29.  This  assumes  that  the  volume  scattering  comes  entirely 
from  the  underlying  sand  sediment.  The  reasoning  behind  this  modeling  comparison 
follows  from  how  the  inversion  is  handling  the  mud  layer.  The  inversion  in  its  current 
state  of  development  treats  the  seafloor  as  a  half-space  and  is  not  currently  equipped  to 
deal  with  a  layered  medium.  As  a  result,  the  processing  interprets  the  first  acoustic  return 
as  coming  from  the  sand/mud  interface,  not  the  mud/water  interface,  and,  in  the  first 
stage  of  the  processing,  the  volume  scattering  is  associated  with  this  sand  layer. 

In  this  case,  the  scattering  strength  determined  from  the  density  fluctuations  is  under 
predicting  the  volume  scattering  by  roughly  10  dB  across  the  entire  band  and  it  is 


f  =  300  kHz 


f  =  400  kHz 


Range  (m) 


Figure  33:  Results  of  a  modified  inversion  which  treats  the  top  15  cm  of  the  half¬ 
space  as  a  layer  with  different  point  volume  scatterers  from  those  scatterers  helow 

the  layer. 
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difficult  to  say  whether  the  inversion  is  failing  or  if  the  seattering  model  used  to  represent 
ground  truth  is  ineomplete.  As  opposed  to  a  sand/water  interfaee  where  the  boundary  is 
elearly  defined,  the  mud  to  sand  transition  is  mueh  smoother  as  can  be  seen  in  Figure  28, 
and  neither  the  seattering  predietion  nor  the  inversion  eurrently  aeeounts  for  this.  This 
interfaee  is  also  very  rough,  5-10  dB  rougher  than  TREX13,  and  it  beeomes  difficult  to 
discriminate  between  interfaee  scattering  and  volume  scattering  along  this  interface  sinee 
the  boundary  is  no  longer  elear.  The  inversion  may  be  treating  this  interfaee  as  volume 
seattering  where  the  seattering  based  on  the  density  fiuetuations  is  not  taking  this  into 
aeeount.  This  hypothesis  is  also  eonsistent  with  the  speetral  inversions  shown  in  Figure 
35  where  the  inversion  slightly  under  estimates  the  sand/mud  interfaee  roughness 
speetrum. 

This  hypothesis  is  strongly  supported  by  inversions  that  we  performed  using  a  more 
eomplex  model.  This  model  estimates  volume  scattering  strengths  for  two  regions,  a  15- 
em  layer  starting  at  the  mud-sand  interface  and  the  half  spaee  below.  This  half  spaee 
eontains  the  regions  in  whieh  the  heterogeneity  measurements  leading  to  the  red  "ground 
truth"  values  in  Figure  31  and  Figure  32  were  made.  For  the  more  eomplieated 
inversions,  volume  scattering  strengths  for  the  15-em  layer  are  similar  to  those  obtained 
from  the  simpler  inversion,  but  those  for  the  lower  half-spaee  agree  with  ground-truth 
within  measurement  uneertainties  (Figure  33).  Thus,  the  inversion  agrees  with  ground 
truth  when  it  isolates  the  region  in  whieh  the  ground-truth  measurements  were  made. 

At  this  stage  it  is  clear  that  the  more  eomplieated  inversion  provides  more  realistie 
results,  but  it  is  not  elear  that  the  resulting  eomplieation  is  an  advanee  in  praetieal  terms. 
First,  the  more  eomplieated  inversion  used  a  "man  in  the  loop",  setting  the  15-em  layer 
thiekness.  In  a  praetieal  implementation,  this  setting  eould  likely  be  automated,  using  the 
eeho-sounder  mode  deseribed  earlier.  Seeond,  as  noted  earlier,  the  simpler  inversion  is 
eonsistent  with  the  physies  used  in  typieal  applieations,  sueh  as  PC  SWAT,  and  will  give 
the  eorreet  reverberation  levels,  even  if  some  of  the  inverted  parameters  are  somewhat 
non-physieal.  Our  experienee  indicates  that  this  situation  will  only  arise  in  difficult 
environments  like  that  of  BayExl4,  and  not  in  the  majority  of  cases. 
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The  results  from  the  lOI  shown  in  Figure  34  and  Figure  36  indieate  that  the  sediment  is 
mud;  The  density  ratio,  porosity,  and  sound  speed  ratio  all  compare  very  well  with  the 
values  measured  in  the  mud  layer.  The  loss  parameter  in  Figure  34  does  not  show  the 
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Figure  34:  Loss  Parameter  determined  from  stage  1  of  the  inversion,  from  the 
inversion  lOI  result,  and  the  attenuation  array. 
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Figure  35:  Roughness  Spectra  along  the  BayExl4  track  compared  to  the  spectra 
measured  using  the  SLS,  Each  panel  corresponds  to  a  different  wavenumber. 
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same  agreement  but  the  attenuation  measurements  in  both  the  sand  and  mud  are 
unreliable  and  the  values  given  from  the  lOI  regressions  are  consistent  values  for  other 
mud  sites.  Again  this  seems  surprising  given  that  the  inversion  is  seeing  the  sand/mud 
layer  as  the  start  of  the  half  space,  but  the  conductivity  probe  measurements  indicate  that 
this  boundary  is  not  abrupt  but  that  there  may  be  a  transition  layer  where  the  sand-to-mud 
ratio  gradually  increases  across  the  boundary.  This  kind  of  transition  layer  has  been 
shown  to  greatly  reduce  the  reflection  coefficient  [23],  particularly  near  normal 
incidence,  which  could  lead  to  a  reduction  in  the  lOI. 

It  appears  that  both  the  inversion  and  the  scattering  models  fail  to  account  for  the 


Range  (m) 


Figure  36:  Sound  Speed  Ratio  (upper),  Porosity  (middle),  and  Density  Ratio  (lower) 
for  BayExl4  determined  using  the  lOI  determined  from  the  inversion.  The 
measured  values  are  also  shown. 
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complexity  of  the  environment.  The  inversion  treats  this  layered  medium  as  a  half-spaee, 
with  no  layer,  and,  as  a  result,  is  ignoring  the  mud  layer.  Treating  the  medium  as  a  half¬ 
space,  it  assumes  that  the  sand/mud  boundary  is  abrupt  and  also  ignores  the  transition 
layer  that  is  clearly  visible  in  the  conduetivity  probe  data.  This  leads  to  low  values  for 
the  reflection  coefficient  and  hence  the  lOI  produced  at  this  transition.  While  it  fails  to 
aceount  for  this  smooth  transition,  it  does  seem  to  treat  the  variation  in  depth  of  the 
sand/mud  interfaee  as  volume  scattering. 

This  complexity  could  be  captured  by  the  inversion  by  introdueing  more  free  parameters 
into  the  first  stage  of  proeessing  as  seen  in  Figure  33.  The  mud  layer  could  be  accounted 
for  by  allowing  the  possibility  of  a  layer  by  introducing  depth-dependent  volume 
scattering  and/or  additional  2D  distributions  of  point  interfaee  scatterers  at  arbitrary 
depths.  This  would  greatly  inerease  the  number  of  parameters  that  the  inversion  would 
need  to  seareh  over  both  inereasing  the  eomputation  time  and  the  uneertainties  in  the 
output.  Ignoring  these  diffieulties  for  the  moment,  suppose  that  the  inversion  works  and 
it  eorreetly  identifies  aspeets  of  the  sediment  sueh  as  the  mud  layer  and  the  transition 
layer  from  the  mud  into  the  sand.  For  the  applications  that  are  of  interest  to  SERDP, 
these  outputs  would  need  to  be  used  as  inputs  to  models  such  as  PC  SWAT  or  into  target 
seattering  models.  These  models,  however,  are  not  currently  designed  to  handle  this  level 
of  complexity  and  as  result  would  fail  to  eorreetly  capture  the  aeoustie  response  of  the 
environment.  It  is  reasonable  to  eonelude  that  the  complexity  of  the  inversion  algorithm 
should  reflect  the  complexity  of  the  models  that  rely  on  its  output. 

If  the  simplicity  of  the  inversion  algorithm  is  maintained  even  for  this  complex 
environment,  are  the  outputs  of  the  inversion  correct?  The  inversion  fits  the  MBES  time 
series  so  the  outputs  are  those  required  to  produce  the  best  fit  of  the  data  for  the 
frequeney  band  and  the  ineident  angles  used  in  the  inversion.  As  a  result,  using  these 
outputs  with  a  eode  sueh  as  PC  SWAT  should  produee  the  eorreet  reverberation  levels 
under  the  same  conditions.  The  question  remains  as  to  whether  the  results  from  PC 
SWAT  would  still  be  eorreet  if  it  were  to  model  smaller  grazing  angles  or  frequencies 
using  the  inversion  output.  To  properly  address  this  question  would  require  that  the 
models  used  to  evaluate  the  inversion  performance  to  account  for  the  complexity  of  the 
environment  and  produce  the  effective  response  of  the  sediment.  While  there  are  tools 
available  to  extend  our  seattering  models  to  capture  this  eomplexity,  due  to  the  difficulty 
of  the  problem  and  amount  of  time  this  effort  would  require,  it  could  not  be  accomplished 
as  a  part  of  this  projeet. 

Anecdotally,  there  is  evidence  that  the  BayExl4  inversion  did  produee  parameters  that 
eorreetly  capture  the  aeoustie  response  of  the  medium  at  lower  grazing  angles  and 
frequencies.  In  performing  data/model  comparisons  for  the  target  scattering  experiments 
that  took  plaee  during  BayExI4,  Dr.  Steve  Kargl  and  Dr.  Aubrey  Espana  have  found  the 
best  eomparisons  with  the  data  occur  when  the  sediment  density  is  assumed  to  be  1.3 
g/em  eonsistent  with  the  results  of  the  inversion  [24].  This  data  is  for  scattering  from 
targets  that  are  resting  on  the  sand/mud  interface  and  for  acoustic  signals  incident  at 
shallow  grazing  angles  and  frequencies  less  than  30  kHz.  The  model  in  this  ease  assumes 
a  half-spaee  and  neglects  the  mud  layer. 
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Conclusions  and  Implications  for  Future  Research/Implementation 

The  goal  of  this  project  was  to  develop  a  high  frequency,  physics-based  sonar  inversion 
and  test  its  performance  in  a  series  of  field  experiments.  The  first  of  these  experiments, 
GulfExll,  which  took  place  prior  to  the  start  of  this  project,  was  hampered  by  the 
presence  of  fish,  problems  with  the  environmental  measurements,  and  issues  with 
saturation  of  the  7125  returns.  Despite  these  issues,  the  data  provided  a  starting  point  for 
the  inversion  development  and  the  opportunity  to  learn  what  was  needed  and  what  should 
be  avoided  in  the  next  two  experiments.  These  lessons  were  successfully  applied  in 
TREX13  and  BayExl4,  and  data  sets  were  collected  of  sufficient  quality  to  perform  the 
inversion  evaluations  that  were  needed. 

Eor  TREX13,  the  MBES  had  to  be  mounted  on  a  pole  alongside  the  ship,  which  led  to 
instabilities  and  artifacts  in  some  of  the  data.  As  a  result,  we  had  to  develop  a  rejection 
criterion  to  remove  these  spurious  data  points,  thus  reducing  our  data  coverage  of  the  site. 
This  turned  out  not  to  be  an  issue  for  a  majority  of  the  data  products  from  the  inversion 
since  they  are  frequency  independent  and  can  be  averaged  across  the  entire  frequency 
band.  The  exception  was  the  volume  scattering  strength,  which  is  frequency  dependent, 
and  the  results  were  sparse  both  along  track  and  in  frequency.  Placing  the  sonar  in  the 
moon  well  of  the  vessel  during  BayExl4  significantly  improved  the  quality  and  stability 
of  the  data  and  we  would  expect  a  similar  improvement  in  data  quality  for  the  TREX13 
site  if  the  sonar  had  been  placed  in  the  well. 

At  the  TREX13  experiment  site,  there  were  two  seabed  regions  covered  by  the  MBES 
track;  A  region  which  was  composed  sand  with  a  shell  and  shell  hash  fraction  that 
increased  with  range  and  a  region  beyond  200  m  along  the  track  that  was  composed  of 
mud  and  sand  in  unknown  proportions  up  until  350  m  where  it  changed  back  to  a  shell 
and  sand  mix.  The  goal  of  this  field  test  was  to  test  the  inversion  using  a  sand 
environment  and  hence  the  first  200  m  were  the  focus  of  the  environmental 
characterization.  In  this  region,  the  roughness  power  spectrum  obtained  by  inversion 
compared  well  with  the  ground  truth  measurements  and  inversion  slightly  overestimated 
the  values  for  most  of  the  lOI  derived  quantities.  This  overestimation  is  consistent  with 
an  increase  in  the  reflection  coefficient  due  to  the  large  quantity  of  shell  hash  lying  on  the 
seafloor  in  this  region. 

The  loss  parameter  was  determined  during  both  the  first  stage  of  the  inversion  and  from 
the  lOI  found  in  the  second  stage  of  the  inversion.  The  loss  parameter  from  the  second 
stage  compared  well  with  the  ground  truth  measurements  while  the  loss  parameter  from 
the  first  stage  was  a  significantly  lower.  This  indicates  that  the  inverted  loss  accounts  for 
only  the  intrinsic  attenuation,  while  the  lOI  loss  captures  both  the  intrinsic  and  scattering 
losses.  The  volume  scattering  strengths  determined  from  the  inversion  underestimate  the 
model  predictions,  but  this  may  be  a  failure  of  the  model  to  account  for  the  significant 
levels  of  multiple  scattering  that  should  be  expected  to  occur  in  this  environment.  There 
are  also  fluctuations  in  the  inverted  volume  scattering  strength  along  the  track  that  may 
be  due  to  variations  in  the  environment  or  instabilities  in  the  sonar  mounting.  The 
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relative  importance  of  these  two  sources  of  variation  cannot  be  determined  using  this  data 
set. 

The  BayExl4  environment  was  arguably  an  order  of  magnitude  more  complicated  than 
the  TREX13  site  in  that  there  was  a  mud  layer  over  a  mixed  sand  mud  half-space  and  a 
transition  layer  that  had  significant  depth  fluctuations  (characterized  as  roughness  in 
Figure  28).  At  its  current  level  of  sophistication,  the  inversion  treated  this  environment 
as  a  half-space  and  as  a  result  the  outputs  were  effective  values  based  on  this 
approximation.  These  effective  values  did  not  necessarily  reflect  the  ground  truth 
measurements,  but  there  is  reason  to  believe  that  these  may  be  the  correct  parameters  to 
use  in  a  model  such  as  PC  SWAT.  Due  to  the  complexity  of  the  environment,  at  this  time 
we  can  only  hypothesize  what  mechanisms,  such  as  the  transition  layer  and  volume 
scattering  at  the  sand/mud  interface,  need  to  be  included  in  the  scattering  models  to 
perform  a  proper  evaluation  of  the  inversion's  performance. 

Despite  this  hurdle  to  fully  evaluating  the  inversion  performance  for  BayExl4,  the 
inversion  has  reached  a  level  of  maturity  with  this  field  test  that  puts  it  at  the  cutting  edge 
of  remote  seabed  characterization  particularly  for  those  applications  that  are  of  most 
interest  to  the  UXO  detection  problem.  Taking  a  physics-based  approach  insures  that  the 
output  of  the  inversion  addresses  the  needs  of  modeling  codes  and  simulations  that 
incorporate  similar  physical  models.  Any  inversion  becomes  unstable  as  the  complexity 
of  its  underlying  model  is  increased.  This  complexity  can  be  constrained  by  the 
maintaining  the  same  level  of  complexity  of  the  scattering  and  simulation  codes  that  the 
inversion  is  being  designed  to  support.  Elsing  the  half-space  model,  the  inversion  was 
able  to  be  applied  seamlessly  to  both  the  TREX13  and  the  BayExl4  MBES  data, 
producing  stable  results  in  these  complex  environments  while  still  providing  meaningful 
outputs. 

Focusing  on  the  application  of  the  inversion  data  products  also  suggests  a  more 
meaningful  approach  to  evaluating  the  inversion  performance: 

1 .  Collect  MBES  data  and  invert  the  data. 

2.  Collect  a  second  set  of  acoustic  data  in  the  same  area  using  a  different  sonar  (side 
scan,  SAS,  etc.). 

3.  Simulate  the  second  set  of  acoustic  data  using  the  outputs  from  the  inversion  and 
compare  that  simulation  to  the  measured  data. 

This  approach  removes  the  requirement  for  ground  truth  measurements  and  directly 
addresses  the  SERDP  needs.  A  similar  approach  could  be  taken  to  support  modeling  of 
UXO  burial  but  it  is  unclear  to  the  author  at  this  time  what  inputs  are  needed  to  those 
models  and  to  what  level  of  sophistication  those  models  have  evolved. 

Our  experiences  with  data  collection  in  these  three  experiments  also  indicate  that  a  new 
data  collection  approach  is  needed.  In  TREX13  and  BayExl4,  data  was  collected  along 
the  measurement  track  by  driving  the  boat  back  and  forth  and  changing  the  frequency  for 
each  pass.  This  is  time-consuming  and  impractical  for  a  wide  area  survey.  A  better 
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approach  would  be  to  change  frequencies  for  each  ping  during  a  pass  along  a  single 
survey  track.  For  those  data  products  that  are  frequency  independent,  the  data  density 
would  not  reduce.  For  the  volume  scattering  strength,  the  reduction  in  data  density  at 
each  frequency  may  be  offset  by  the  improvements  in  data  quality  that  would  come  from 
collecting  all  of  the  data  in  a  single  pass.  Discussions  with  engineers  at  Teledyne- 
RESON  indicate  that  this  approach  could  be  implemented  with  the  addition  of  an  external 
PC  control  that  would  avoid  modifying  the  existing  suite  of  software  that  runs  the  sonar. 
This  not  only  simplifies  the  implementation  but  also  provides  a  transition  path  where  this 
capability  could  be  added  to  existing  7125s.  Development  of  this  system  would  require 
that  we  move  away  from  the  approach  taken  in  this  project  where  we  have  borrowed 
sonars  from  the  Teledyne -RESON  and  instead  purchase  a  sonar  that  could  be  used 
exclusively  in  developing  and  testing  the  frequency-hopping  approach. 
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